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We calculate the most interesting rare and CP-violating K and B decays in the 
Littlest Higgs model with T-parity. We give a collection of Feynman rules includ- 
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ing v 2 / f 2 contributions that are presented here for the first time and could turn 
out to be useful also for applications outside flavour physics. We adopt a model- 
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independent parameterization of rare decays in terms of the gauge independent 
functions X{, Yi, Z{ (i = K, d, s), which is in particular useful for the study of the 
breaking of the universality between K, and B s systems through non-MFV in- 
teractions. Performing the calculation in the unitary and 't Hooft-Feynman gauge, 
we find that the final result contains a divergence which signals some sensitivity 
to the ultraviolet completion of the theory. Including an estimate of this contribu- 
tion, we calculate the branching ratios for the decays K + — > ir + vv, Kl — > iftvv, 
B s ,d — > B — > X s ^vv, Kl — > tt°£ + £~ and B — > X Sj d^ + £~ , paying particular 

attention to non-MFV contributions present in the model. 

The main feature of mirror fermion effects is the possibility of large modifi- 
cations in rare K decay branching ratios and in those B decay observables, like 
S^a, and -AgL, that are very small in the SM. Imposing all available constraints 
we find that the decay rates for B s ^ — > fi + fj,~ and B — > X S ( jyi> can be enhanced 
by at most 50% and 35% relative to the SM values, while Br(K + — > ^vv) and 
Br(Ki — > ifivv) can be both as high as 5 • 1CP 10 . Significant enhancements of the 
decay rates Kl — > tt°£ + £~ are also possible. Simultaneously, the CP-asymmetries 
S^tf, and A| L can be enhanced by an order of magnitude, while the electroweak 
penguin effects in B — > irK turn out to be small, in agreement with the recent 
data. 



Note added 



An additional contribution to the Z penguin in the Littlest Higgs model with T-parity 
has been pointed out in [1,2], which has been overlooked in the present analysis. This 
contribution leads to the cancellation of the left-over quadratic divergence in the cal- 
culation of some rare decay amplitudes. Instead of presenting separate errata to the 
present work and our papers [3-6] partially affected by this omission, we have presented 
a corrected and updated analysis of flavour changing neutral current processes in the 
Littlest Higgs model with T-parity in [7]. 



1 Introduction 



Rare and CP-violating K and B meson decays will hopefully provide in the coming 
years a new insight into the origin of the hierarchy of quark masses and their hierarchical 
flavour and CP- violating interactions. While the presently available data on these decays 
give a strong indication that the CKM matrix [8] and more generally minimal flavour 
violation (MFV) [9-11], encoded entirely in Yukawa couplings of quarks and leptons, is 
likely to be the dominant source of flavour and CP violation, there is clearly still room 
left for contributions governed by non-MFV interactions, in particular those with new 
CP-violating phases. 

A prime example of a model with non-MFV interactions is a general MSSM in which 
non-MFV contributions originate in squark mass matrices that are not aligned with the 
quark mass matrices. Extensive analyses of the impact of such non-MFV contributions 
on particle-antiparticle mixing and rare K and B decays in the MSSM with and without 
i?-parity have been presented in the literature. The basic strategy proposed in [12] is to 
constrain the non-diagonal entries (Sij) ab of the squark mass matrices given in the so- 
called super-CKM basis, in which all neutral gauge interactions and the quark and lepton 
mass matrices are flavour diagonal. These studies are complicated by the fact that in 
addition to (Sij) ab also new fermion and scalar particle masses appear as free parameters. 
The situation in this respect may improve significantly in the coming years provided the 
supersymmetric particles will be discovered at LHC and their masses measured at LHC 
and later at ILC. 

While supersymmetry [13] appears at present to be the leading candidate for new 
physics beyond the Standard Model (SM), Little Higgs models [14,15] appear as an 
interesting alternative. Also here, new particles below 1 TeV, or slightly above, with sig- 
nificant contributions to electroweak precision studies and FCNC processes are present. 
Among the most popular Little Higgs models is the so-called Littlest Higgs model [16]. 
In this model in addition to the SM particles, new charged heavy vector bosons (W^), 
a neutral heavy vector boson (Z H ), a heavy photon (Ah), a heavy top quark (T + ) and 
a triplet of scalar heavy particles (<&) are present. 

In the original Littlest Higgs model (LH) [16], the custodial SU(2) symmetry, of 
fundamental importance for electroweak precision studies, is unfortunately broken al- 
ready at tree level, implying that the relevant scale of new physics, /, must be at least 
2 — 3 TeV in order to be consistent with electroweak precision data [17-23]. As a conse- 
quence, the contributions of the new particles to FCNC processes turn out to be at most 
10 — 20% [24-28], which will not be easy to distinguish from the SM due to experimental 
and theoretical uncertainties. In particular, a detailed analysis of particle-antiparticle 
mixing in the LH model has been published in [24] and the corresponding analysis of 
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rare K and B decays has recently been presented in [28]. 

More promising and more interesting from the point of view of FCNC processes 
is the Littlest Higgs model with a discrete symmetry (T-parity) [29] under which all 
new particles listed above, except T+, are odd and do not contribute to processes with 
external SM quarks (T-even) at tree level. As a consequence, the new physics scale / 
can be lowered down to 1 TeV and even below it, without violating electroweak precision 
constraints [30]. 

A consistent and phenomenologically viable Littlest Higgs model with T-parity (LHT) 
requires the introduction of three doublets of "mirror quarks" and three doublets of 
"mirror leptons" which are odd under T-parity, transform vectorially under SU (2) L and 
can be given a large mass. Moreover, there is an additional heavy T_ quark that is odd 
under T-parity [31]. 

In the first phenomenological study of FCNC processes in the LHT model by Hubisz 
et al. [32], the impact of mirror fermions on the neutral meson mixing in the K, B and D 
systems has been studied in some detail. As described in that paper, in the LHT model 
there are new flavour violating interactions in the mirror quark sector, parameterized by 
two CKM-like mixing matrices Vna and Vh u , relevant for processes with external light 
down-type quarks and up-type quarks, respectively. It turns out that the spectrum of 
mirror quarks must be generally quasi-degenerate, but there exist regions of parameter 
space, where only a loose degeneracy is necessary in order to satisfy constraints coming 
from particle-antiparticle mixing. 

In [33] we have confirmed the analytic expressions for the effective Hamiltonians for 
K° — K°, B® — B® and B® — B° s mixings presented in [32] and we have generalized 
the analysis of these authors to other quantities that allowed a deeper insight into the 
flavour structure of the LHT model. While the authors of [32] analyzed only the mass 
differences AMk, AM^, AM,, AMo and the CP-violating parameter ek, we included in 
our analysis also other theoretically cleaner quantities: the CP asymmetries Acp(Bj — > 
ipKs), Aqj>(B s — > ip<f)) and A^l, and the width differences Ar^ s . We have also calculated 
Br(B — > X s ^,fi) and the corresponding CP asymmetries. Moreover we emphasized that 
while the original LH model belongs to the class of models with MFV, this is certainly 
not the case of the LHT model where the presence of the matrices Vnd and Vhu in the 
mirror quark sector introduces new flavour and CP-violating interactions that could have 
a very different pattern from the ones present in the SM. 

The main messages of our analysis in [33] are as follows: 

• The analysis of the mixing induced CP asymmetries Ao^iB^ — > ipK$) and Aqp{B s 
— > ip<f)) illustrates very clearly that with mirror fermions at work these asymme- 
tries do not measure the phases — j3 and — j3 s of the CKM elements Vu and V ts , 
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respectively. 

• This has two interesting consequences: first, the possible "discrepancy" between 
the values of sin 2(5 following directly from A C p(B d — > i>K s ) and indirectly from 
the usual analysis of the unitarity triangle involving AM^ S , ek and IV^/K&I can 
be avoided within the LHT model. Second, the asymmetry Acp{B s — > ijxf)) can be 
significantly enhanced over the SM expectation ~ 0.04, so that values as high as 
0.3 are possible. Moreover, the asymmetry A$ L can be enhanced by an order of 
magnitude. 

• We also found that the usual relation between AM d /AM s and |\4i/V^| character- 
istic for all models with MFV is no longer satisfied. 

• Calculating BriB — > X s '-f) in the LHT model for the first time we have found 
that in spite of large effects in AF = 2 processes considered in [32,33], it can be 
modified by at most ±4% relative to the SM value [34, 35] . This is welcome as the 
SM branching ratio agrees well with the data. Also new physics effects in B — > Xjy 
and in the CP asymmetries in both decays are at the level of a few percent of the 
SM values. 

The beauty of the LHT model, when compared with other models with non-minimal 
flavour violating interactions, like general MSSM, is a relatively small number of new 
parameters and the fact that the local operators involved are the same as in the SM. 
Therefore the non-perturbative uncertainties, present in certain quantities already in 
the SM, are the same in the LHT model. Consequently the departures from the SM 
are entirely due to short distance physics that can be calculated within perturbation 
theory. In stating this we are aware of the fact that we deal here with an effective field 
theory whose ultraviolet completion has not been specified, with the consequence that 
at a certain level of accuracy one has to worry about the effects coming from the cut-off 
scale A ~ 47r/. 

As pointed out recently in [28], such effects are signaled by left-over logarithmic 
divergences in the final result for FCNC amplitudes and they appear as poles 1/e when 
dimensional regularization is used. It turns out that such divergences are absent in 
particle-antiparticle mixing and B — > X s>l f) decays both in the LH and the LHT model. 
On the other hand, they are present in Z°-penguin diagrams in the LH model considered 
in [28] . As we will see in the present paper, the imposition of T-parity eliminates such 
divergences from the T-even sector as the diagrams containing these divergences are 
forbidden by T-parity. However, we find that the mirror fermion contributions to AF = 1 
processes contain such logarithmic divergences and their effects on rare decays has to be 
taken into account following the steps of [28]. 



3 



In the present paper we extend our analysis of [33] to include all prominent rare K 
and B decays. In particular we calculate the LHT contributions from both T-even and 
T-odd sectors to the one- loop functions X, Y and Z [36] that encode the short distance 
contributions to the rare decays in question. The new properties of these functions 
relative to MFV are: 

• They are complex quantities, as in the phenomenological analysis in [37], while in 
contrast to the real X, Y and Z of MFV. 

• They are different for K, B d and B s systems in contrast to [37], where these 
functions were the same for the three systems in question. Thus in the present 
model we deal really with nine short distance functions. But as we calculate them 
in a specific model, interesting correlations between observables in rare K, B d and 
B s decays are still present, although their structure differs from the correlations in 
MFV models [38-40] and in the phenomenological non-MFV analysis in [37]. 

Our paper is organized as follows. In Section 2 we review those ingredients of the LHT 
model that are of relevance for our analysis. In particular we summarize in an appendix 
the relevant Feynman rules that go beyond those presented in [32,41,42]. Section 3 is 
devoted to the master formulae for rare decays in models with non-MFV interactions. 
This presentation goes beyond the LHT model and the formulae given in this section are 
general enough to be used in other models that contain non-MFV contributions. This 
section is basically a generalization of the formulae in [37] to include the breakdown 
of universality in the functions X, Y and Z as mentioned above. In Section 4 we 
calculate the new contributions to X, Y and Z coming from the T-even sector. The 
results from the analysis of the LH model in [28] were very helpful in obtaining these 
results. In Section 5, the most important section of our paper, we calculate the T-odd 
contributions to X and Y in the 't Hooft-Feynman and unitary gauge, obtaining a gauge 
invariant result for these functions. Similarly to [28] we identify logarithmic divergences 
that are gauge independent. We recall the interpretation of these divergences and we 
analyze their impact on the functions X and Y in Section 6. In Section 7 we calculate 
the function Z that contributes to the B — > X s i + i~ decay. In Section 8 we study the 
decay K L — > ir°l + l~. The benchmark scenarios for our numerical analysis are described 
in Section 9. In Section 10, the numerical analysis of the relevant branching ratios 
including T-even and T-odd contributions is presented in detail. Of particular interest 
is the pattern of the violations of several MFV relations [38-40] that can be tested 
experimentally one day. This analysis takes into account the constraints found from the 
analysis of AF = 2 processes and B — » X s 7 presented in [33] and from Br(B — > X s £ + l~) 
calculated here. It can be considered as the first global analysis of FCNC processes in 
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the LHT model. In this section we also comment on B — > txK decays. In Section 11 we 
conclude our paper with a list of messages resulting from our analysis and with a brief 
outlook. Few technical details and the Feynman rules can be found in the appendices. 



2 General Structure of the LHT Model 

2.1 Gauge and Scalar Sector 

The original Littlest Higgs model [16] is based on a non-linear sigma model describing the 
spontaneous breaking of a global SU (5) down to a global SO (5). This symmetry breaking 
takes place at the scale / ~ O(TeV) and originates from the vacuum expectation value 
(VEV) of an SU(5) symmetric tensor E, given by 
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(2.1) 



The generators T a of the unbroken SO (5) symmetry fulfill T a S + S (^ a ) T = 0, 
whereas the broken generators X a of SU(5)/SO(5) satisfy X a S - S (X a ) T = 0. The 
symmetry breaking can thus be considered as a Z 2 automorphism under which T a i-)- T a 
and X a i->- — X a . This automorphism is the basic motivation for the way T-parity is 
implemented, as discussed later. 

The mechanism for this symmetry breaking is unspecified, therefore the Littlest Higgs 
model is an effective theory, valid up to a scale A ~ 47r/. From the SU{h) / SO{h) 
breaking, there arise 14 Nambu-Goldstone bosons \ a which are described by the "pion" 
matrix II, given explicitly by 
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Here, H = (-in+fy/Q,, (v + h + ivr°)/2) T plays the role of the SM Higgs doublet, i.e. 
h is the usual Higgs field, v = 246 GeV the Higgs VEV, and tt ± , it are the Goldstone 
bosons associated with the spontaneous symmetry breaking SU(2) L x U(1) Y — > U(l) em . 
The fields rj and 00 are additional Goldstone bosons, and $ is a physical scalar triplet, 
as described below. 
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diag(3,3,-2,-2,-2), 
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diag(2,2,2,-3,-3), 



(2.4) 



The low energy dynamics of the symmetric tensor E is then described by 

E = e m / f Z e mT / f = e^Eo . (2.3) 

An [SU(2) x x [SU(2) x U(l)} 2 subgroup of the SU(5) global symmetry is 

gauged, with the generators 

(a a 0\ 

, Yi = 
\0 0/ 

/o \ 

, Y 2 

\o o -d a 7 

and the corresponding gauge bosons W^, B^, W 2 ^, B%. Having the Z 2 automorphism 
T a I—)- T a and X a h-> — X a in mind, one implements T-parity by assigning a symmetry 
transformation to all fields, such that the theory is symmetric under T-parity. A natural 
way to define the action of T-parity on the gauge fields is 

W? ^ W% , BiO B 2 . (2.5) 

In other words, the action of T-parity exchanges the two SU(2) x U(l) factors. An 
immediate consequence of this definition is that the gauge couplings of the two SU (2) x 
£7(1) factors have to be equal. 

The gauge boson T-parity eigenstates are given by 
W? + W£ „ B 1 + B 2 
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(2.7) 



V2 ' V2 
where "L" and U H" denote "light" and "heavy", respectively. 

Now, the VEV E breaks the gauge symmetry to the diagonal T-even SU{2) x U(l), 
which is identified with the SM gauge group. The breaking SU(2)l x U(1)y — > U(l) em 
then takes place via the usual Higgs mechanism. Finally, the mass eigenstates are given 
at 0(vyp) by 
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Z L = cos6wW L — syhOwBl , 
At = sin 6 w Wl + cos 6 W B L , 



Z H = W^ + x H —B H , 
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(2.9) 
(2.10) 
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where 9w is the usual weak mixing angle and 

5gg' 

From these breakings, the T-odd set of SU(2) xU(l) gauge bosons acquires masses, 
given at 0(v 2 /f 2 ) by 

Mw B = fg(l-j^, M Zh =M Wh} M AH = ^=(^l-^y (2.12) 

The masses of the T-even gauge bosons are generated only through the second step of 
symmetry breaking. They are given by 

m "> = t O-tif)' ^^^X 1 -^)' Ma ^ (l (2 ' 13) 

Note that T-parity ensures that the custodial relation Mw L = Mz L cos9w is exactly 
satisfied at tree level. 

In order to ensure that $ , u , 77 are odd under T-parity, whereas the SM Higgs doublet 
H is even, one makes the following T-parity assignment: 

Il^-firm, where Q — diag(l, 1, -1, 1, 1) . (2.14) 

As mentioned above, $ is an additional scalar triplet in a symmetric representation 
of SU{2)l. Its mass is given by 



= y/2m H -, (2.15) 



v 

with mu being the mass of the SM Higgs scalar. As pointed out in [30] , mu in the LHT 
model can be significantly higher than in supersymmetry. In Appendix A we show that 
$ has only negligible effects on the decays studied in the present paper. 

The fields u ± 1 u° and r\ are the Goldstone bosons associated with the breaking of 
[SU(2) x U(l)]i x [SU(2) x U(l)] 2 to its diagonal subgroup. They are thus eaten by the 
heavy gauge bosons W^, Z H and A H) respectively. 1 



2.2 Fermion Sector 

A consistent and phenomenologically viable implementation of T-parity in the fermion 
sector requires the introduction of mirror fermions [31]. We embed the fermion doublets 

lr To be exact, <£>, 10 and r\ mix with each other at 0(v 2 / f 2 ) [30] and it is a linear combination of the 
fields that is eaten. See Appendix B for details. 
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into the following incomplete representations of £77(5), and introduce a right-handed 
SO (5) multiplet 
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Under T-parity these fields transform as 

Thus, the T-parity eigenstates of the fermion doublets are given by 
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(2.18) 



(2.19) 



5sm are the left-handed SM fermion doublets (T-even), and qn are the left-handed mirror 
fermion doublets (T-odd). The right-handed mirror fermion doublet is given by ipR. 
The mirror fermions can be given O(f) masses via a mass term 



£ mirror = _ K .. f + ^Eo^ 1 "^) *i , 



(2.20) 



where we sum over the generation indices i,j = 1, 2, 3 and £ = e lU ^ is needed to make 
^mirror SU(5) invariant. 

The mirror fermions thus acquire masses, given by [32] 

„,2 



m 



Hi 



m m 1 - 



(2.21) 
(2.22) 



where the eigenvalues of the mass matrix k. 

The additional fermions ipn and xr can be given large Dirac masses by introducing 
additional fermions, as described in detail in [31,41]. In what follows, we will simply 
assume that they are decoupled from the theory. 



2.3 Yukawa Sector 

In order to cancel the quadratic divergence of the Higgs mass coming from top loops, an 
additional heavy quark T + is introduced, which is even under T-parity and transforms, 
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to leading order mv/f, as a singlet under SU(2)l- The implementation of T-parity then 
requires also a T-odd partner T_, which is an exact singlet under SU (2)i x SU(2) 2 and 
therefore does not contribute to the decays in question (see [33] for details). 
The Yukawa coupling for the top sector is then given by [41,43] 
1 



--top 



2^2^f € ijk e xy (Ql)i{E)jx('E')ky ~ (Q2^o)i(^)jx('E')ky 
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where 
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(2.23) 



(2.24) 



the superscript in denotes the third quark generation, and £ = S J1S'QE is the 
image of £ under T-parity. 
Note that under T-parity 

Qi -X Q 2 , f lfl -t^ , u R ^u 3 R , (2.25) 
so that the T-parity eigenstates are given by 

t'i T t'2 ,i _ t'w T t' 2 R 

(-4- 



t' 



(2.26) 



±_ V2 ' v/2 
As i'_ and t'_ R do not mix with the mirror fermions at tree level, the mass eigenstates, 
denoted by (TL) L and (T_) fi for the left and right-handed part, respectively, are simply 
given by 

(T_) L = f_, (T_)fl = f_ fl . (2.27) 

However, the T-even eigenstates mix with each other, so that the mass eigenstates 
of the top quark t and its heavy partner T + are given by 

th = c L (<?sm)i - s L t' + , (T+) L = s L (q S u)i + c L t' + , (2.28) 

t R = c R u R - s R t' +R , (T+) R = s R u 3 R + c R t' +R , (2.29) 

where (<?sm)i denotes the upper component of the left-handed SM quark doublet, and 
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(2.32) 
(2.33) 
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with 
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(2.34) 



A? + A|' 

This mixing leads to a modification of the top quark couplings relatively to the SM, as 
can be seen in the Feynman rules given in Appendix B. 
The masses are then given by 
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As the Yukawa couplings of the other SM quarks are small, there is no need to 
introduce additional heavy partners to cancel their quadratically divergent contribution 
to the Higgs mass. Thus the Yukawa coupling for the other up-type fermions is simply 
given by 2 



-up 



2\/2 

and their masses are given by 



^uftijktxy {Ql)i{^)jx{^)ky — (Q2So)j(S) Ja .(E) fc j / U R + tl.C, (2.38) 
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For the down-type Yukawa term, we take [43] 

^ 2 ) x {T) iy {T) jz X - {^ 1 H Q ) x {t) iy {t) jz x\ d R + h.c. , 



r - J^l f 



(2.40) 



where we sum over i, j — 1,2 and x,y, z — 3, 4, 5, and X = (E 33 ) 1 / 4 has been introduced 
in order to make £down gauge invariant. Note that here 
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2 Strictly speaking, all three generations have to be included in the top Yukawa term, where Ai then 
becomes the usual 3x3 Yukawa coupling matrix, as discussed in detail in [44]. However, as found there, 
the mixing of T + with u, c quarks is experimentally highly constrained, so we can safely neglect it. For 
simplicity, we thus write the Yukawa term for each generation separately, and later include the CKM 
mixing "by hand" in the Feynman rules given in Appendix B. 
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i. e. no insertion of a 2 is needed. From this Yukawa term, we obtain the down-type 
quark masses to be 



>»'< x > ( 1 Tip 



1,2,3) 



(2.42) 



Lepton masses are generated in a completely analogous way. 



2.4 Weak Mixing in the Mirror Sector 



As discussed in detail in [32,33], one of the important ingredients of the mirror sector is 
the existence of four CKM-like unitary mixing matrices, two for mirror quarks and two 
for mirror leptons: 
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They satisfy 



vLvh, = v, 



Hu v Hd ^CKM j v Hv v Hi ^PMNS i (2.44) 

where in Vpmns [45] the Majorana phases are set to zero as no Majorana mass term has 
been introduced for the right-handed neutrinos. The mirror mixing matrices in (2.43) pa- 
rameterize flavour violating interactions between SM fermions and mirror fermions that 
are mediated by the heavy gauge bosons Wh, Z h and A H . The notation in (2.43) indi- 
cates which of the light fermions of a given electric charge participates in the interaction. 
Feynman rules for these interactions are given in Appendix B. 

In the course of our analysis it will be useful to introduce the following quantities [33] : 



AK) _ irtisin, 
Si — V Hd V Hd J 



tisirid 



: (d) 



i/*ibirid 
V Hd V Hd J 



Hd v Hd 



1,2,3), (2.45) 



that govern K, B d and B s decays, respectively. 

Following [46] we will parameterize Vud generalizing the usual CKM parameteriza- 
tion, as a product of three rotations, and introducing a complex phase in each of them, 
thus obtaining 
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u 23 
" s 23 e ^ 3 



o \ 



'23 c 



/ 



-23 



-13 



A AS 1 } 
'13 e 



sf 3 e-<«\ / 



"-12 
d iti d 



-13 



Performing the product one obtains the expression 



s\ 2 e i,5 i 2 



-12 



0\ 



(2.46) 



V, 



Hd 



( 



r .d „d 
c 12 c 13 



' 12^23 c 



c 



d „d „d 



12°23°13 
d r .a 

12 c 23°13 



„d „d 
c 12 c 23 



sf 2 c d 13 e 

q(1 r%d r%d 

O i o O OQ OioG 



12°23°13 

d n d c d 
12 L 23*13 



i(5f 3 -5f 2 -5 2 d 3 ) 



d -8f 2 ) 



sf 3 e 

-,d „d 
'23 L 13 C 



"13 

-i<5. d 



„d „d 
c 23 c 13 



\ 
/ 

(2.47) 
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As in the case of the CKM matrix the angles 9fj can all be made to lie in the first quadrant 
with < 5f 2 , > Sf 3 < 27T. The matrix Vh u is then determined through Vh u — ^HdVcxM- 

We point out that in [32] and in the first version of [33] Vna was parameterized in 
terms of three mixing angles and only one phase like Vckm, overlooking the presence of 
two additional phases. The presence of three phases was first pointed out in [46]. In 
short, the reason for the appearance of two additional phases relative to the CKM matrix 
is as follows. Vckm and Vnd are both unitary matrices containing three real angles and 
six complex phases. Varying independently the phases of ordinary up- and down-quark 
states allows us to rotate five phases away from Vckm (an over-all phase change of all 
the quark states leaves Vckm invariant). In rotating phases away from Vnd, one can still 
act on only three mirror states, thus obtaining for Vud a parameterization in terms of 
three mixing angles and three phases. 

The six parameters of Vud have to be determined in flavour violating processes. In [33] 
we have outlined briefly this determination in the context of particle-antiparticle mixing. 
Including rare K, B d and B s decays will further help to determine these parameters. 



3 Rare K and B Decays beyond MFV 
3.1 Preliminaries 

Before presenting in Sections 4 and 5 the details of the calculations of rare K and B 
decays in the LHT model in question, it will be useful to have a general look at rare decays 
within models with new flavour and CP-violating interactions but with the same local 
operators of the SM or more generally of constrained MFV (CMFV) models, as defined 
in [9,40]. While the presentation given below is tailored to the subsequent sections, it 
can easily be adapted to any model of this type. 

It should be emphasized that while the formulae given below bear many similarities 
to the ones given in [37], they differ from the latter ones in the following important 
manner. In [37] a simple beyond-MFV scenario of new physics has been considered 
in which new physics affected only the Z°-penguin function C that became a complex 
quantity, but remained universal for K, Bd and B s decays. In this manner several CMFV 
relations involving only CP-conserving quantities remained valid and the main new effects 
were seen in CP-violating quantities like Br(K L — > -k q vu) and the CP-asymmetries in 
B — > X s £ + l~. In particular, the full system of rare K, B d and B s decays considered in 
this section could be described by three complex functions 

X=\X\e idx , Y = \Y\e i0Y , Z=\Z\e idz , (3.1) 

with correlations between these functions resulting from the universality of the Z°- 
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penguin function C = \C\exp(i$c). As a result the CMFV correlations between ob- 
servables in K, B d and B s were only affected in the cases in which $i played a role. In 
the LHT model the structure of new flavour violating interactions is much richer. Let 
us spell it out in explicit terms. 

3.2 Xi, Yi, Zi functions 

In the CMFV models the new physics contributions enter for all practical purposes only 
through the functions X, Y and Z that multiply the CKM factors A^ 

\\ K) = V* V td , \i d) = V* V td , Aj s) = V* V ts , (3.2) 

for K, Bd and B s systems respectively. 

It will be useful to keep this structure in the LHT model and absorb all new physics 
contributions in the functions X iy Y iy Zi with i — K,d,s defined as follows: 

Xi = X SM + X cven H — 777-X° dd = \Xi\ e l6>x , (3.3) 

AT 

Y % = Y SM + Fcvcn + 4r^° dd = 1^1 e<4 ' ( 3 - 4 ) 

Zi = Z SM + Z CVCQ H — 7KZ° dd = \Zi\ e l9z . (3.5) 

K 

Here Xsm, ^sm and Z$m are the SM contributions for which explicit expressions can 
be found in Appendix C. X even , y^ ven and Z even are the contributions from the T-even 
sector, that is the contributions of T + and of t at order v 2 / f 2 necessary to make the GIM 
mechanism [47] work. The latter contributions, similarly to X SM , Y SM and Z SM , are real 
and independent of % = K, d, s. They can be extracted from [28] and will be given in 
Section 4. The main result of the present paper is the calculation of the functions A"? dd , 
Y° dd and Z° dd , that represent the T-odd sector of the LHT model and are obtained from 
penguin and box diagrams with internal mirror fermions. The details of this calculation 
can be found in Section 5. In what follows we will present the most interesting branching 
ratios in terms of Xj and Y { . The CKM elements that we will use are those determined 
from tree level decays. 

3.3 K -KVV 

Generalizing the formulae in [37] we have 

Br(K + ->■ ti + vv) = k + [f 2 A 4 R 2 \X K \ 2 + 2fP c (x)A 2 R t \X K \ cos/3^ + P c (x) 2 ] , (3.6) 
Br(K L -> tv° vis) = K L f 2 A A R 2 \X K \ 2 sin 2 p$ , (3.7) 
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where [48] 
V ts 



cb 



0.98, k + = (5.08 ±0.17) • HT 11 , k l = (2.22 ± 0.07) • 10~ 10 , (3.8) 
A 2 " 



P c (x) 



1 - — ) P c (x) , P c (x) = 0.42 ± 0.05 , 



(3.9) 



with P c (x) including both the NNLO corrections [48] and long distance contributions [49]. 
Finally 

0$ = 0-0,-0$- (3.10) 

The values of A, R b) (3 and S are collected in Table 1 of Section 10. 
Of particular interest is the relation 



sin2(/3 + v?B d ) = sin2/3f , 



(3.11) 



that for (f Bd — 0, 9% — reduces to the MFV relation of [50,51]. A violation of 
this relation would signal the presence of new complex phases and generally non-MFV 
interactions. In this context the ratio 

2 r 

:1V* 

(3.12) 



Br(K L -> tt°uu) 



Br(K L -» ir°vv] 



SM 



X 



K 



x 



SM 



oK 1 2 



sin ((3 - /3 S 



is very useful, as it is very sensitive to 9$ and is theoretically very clean. 
The most recent SM predictions for the branching ratios read [48] 

Br(K + -> tt+z/z/) = (8.0±l.l)-10- u , Br(K L -> k°uu) = (2.9±0.4) • 10' 11 , (3.13) 

to be compared with the present experimental measurements [52, 53] 

Br(K + -> 7i + uu) = (1.47lJJ°) • 10~ 10 , Br(K L -> vrVz/) < 2.1 • 10~ 7 (90%C.L.) . 

(3.14) 

Recent reviews of the K — > -kvv decays can be found in [54,55]. 
3.4 B s , d --> 

Here, we will mainly be interested in the following ratios 



Br(B, 



fi + fi ) 



Br(B s -> //+// )sm 

Br(B d ^fj. + (j.-) 
Br(B d -> //+//- ) S m 

Br(B d ^fj,+ti-) 
Br(B s -> //+//-) 



SM 



1c 



SM 
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(3.15) 
(3.16) 
(3.17) 



where the departure of the last factor from unity signals non-MFV interactions. In 
obtaining these formulae we assume that the CKM parameters have been determined 
in tree level decays independently of new physics so that they cancel in the ratios in 
question. 

In the LHT model 



AM d _ m Bd B Bi F 2 Bd 
m Bs B B Fl 



where 



AM, 



V t 



Id 



v t 



ts 



c 



B d 



c r 



AM, 



(q = d,s). 



(3.18) 



(3.19) 



(AM ? ) SM 

Consequently, using (3.17) and (3.18), the golden relation between Br(B ds — > fi + fi~) 
and AM d /AM s valid in CMFV models [38] gets modified as follows: 

2 



Br(B s -» fi+fj,-) = B Bd t(B 8 ) AM s 
Br{B d -)• (jL+fi-) B Bs r(B d ) AM d 



r, 



C B d 
Cn 



(3.20) 



with r being generally different from unity. 
The most recent SM predictions read [40] 



Br[B s ->• fj+fT) = (3.35 ± 0.32) • 10~ 9 , Br(B d ->■ fj+fT) = (1.03 ± 0.09) • HT 10 , 

(3.21) 

to be compared with the experimental upper bounds from CDF [56] 



Br(B s ->■ fj+fj,-) < 1 • 10~ 7 , Br(B d ->■ < 3 • 10~ 8 . 



(3.22) 



3.5 B — y X s ^ d uu 

We will also study the theoretical clean decay B — > X S:d uP and look at the ratios 



Br(B 



X s vv) 



Br(B ->■ X s uu) SM 

Br{B ->■ X d i/p) 
Sr(S ->■ X d i/p) SM 

Sr(S ->■ X d vv) 
Br(B -)> X s i/p) 



X s 



SM 



SM 



ts 



(3.23) 
(3.24) 
(3.25) 



with \V td \ and |T4 S | obtained from tree level decays. Note that for X d ^ X s the relation 
of the last ratio to [V^/V^l is modified with respect to MFV models. 
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4 Results for the T-even Sector 



The contribution from the T-even sector, denoted by X cvcn and F e ven in (3.3) and (3.4), 
respectively, can be extracted from [28], where the functions X and Y have been calcu- 
lated in the LH model without T-parity. Imposing T-parity implies 



s = c= s = c = 



V2 



(4.1) 



and vanishing of the diagrams in the classes 1, 2, 4 and 6 in [28]. Moreover there are 
no corrections from the breakdown of custodial symmetry, and the left-over divergence, 
discussed in detail in [28], is also absent. 
We find then 



x, 



x, 



V 

•p 

7 2 



U 3 (x t ,x T ) 
V 3 (x t ,x T ) 



XL Xt 

\-x L % 

XL x t 

l-x L ~& 



(4.2) 
(4.3) 



with the two terms on the r.h.s. coming from class 3 and 5 in [28], respectively. The 
functions Us(x t ,XT) and Vz{x t ,XT) are given in Appendix C. 



5 Results for the T-odd Sector 

5.1 Preliminaries 

In Figs. 1 and 2 we show the diagrams contributing to X° dd in (3.3). Similar diagrams 
but with external charged leptons contribute to Y° dd in (3.4). They can be divided into 
two classes. The first class involves in the unitary gauge and the mirror fermions u l H 
exchanges, while the second class involves Zh, Ah and d H exchanges. In the renormal- 
izable gauges also diagrams with Goldstone bosons have to be included. We will first 
discuss the calculation in the unitary gauge. Subsequently we will turn to the calculation 
in the 't Hooft-Feynman gauge, verifying our result in the unitary gauge. 

5.2 Z^-Penguin Diagrams in the Unitary Gauge 

Only Z L -penguin diagrams contribute to the decays considered because the couplings of 
Z H and A H to vv and vanish due to T-parity. We note that the diagrams with 

internal are fully analogous to the corresponding SM diagrams with internal W^. 
Moreover the diagrams with triple gauge boson vertices vanish in the case of internal 
Ah and Zh contributions. 



16 




Figure 1: Z L -penguin diagrams contributing in the T-odd sector. 



There are two additional features with respect to the SM calculation and the box 
diagram calculation presented in [33] and below: 

• The diagrams in Fig. 1 appear first to be 0(1), that is they are not suppressed by 
v 2 /P. 

• The couplings of mirror fermions to are vectorial (7^) in contrast to the SM 
couplings that have both 7^ and 7^75 components. 

Clearly the 0(1) contributions have to vanish as otherwise it would not be possible to 
decouple the mirror fermions in the limit / — » 00. This is assured by the vectorial 
coupling of Zl to the mirror fermions. The missing of diagrams with triple gauge bo- 
son vertices in the neutral gauge boson case is compensated by the difference between 
d l H Z£d H and u l H Z£u H couplings, so that the charged (W^) and neutral (Zh, Ah) gauge 
boson contributions of 0(1) to the Z^-penguin vanish independently of each other in the 
unitary gauge. 

As the inclusion of v 2 /f 2 corrections to the neutral gauge boson interactions leads 
only to an overall factor multiplying the Zh and Ah contributions, which vanish inde- 
pendently of each other, we find that there is no contribution from mirror fermions to 
■Z^-penguin diagrams in the unitary gauge. The inclusion of v 2 /f 2 corrections to the 
relations between the masses of u H and d l H and to the gauge boson masses does not 
change this result. 
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Figure 2: Box diagrams in the unitary gauge. 

5.3 Box Diagrams in the Unitary Gauge 

In order to simplify the formulae, we will present the results for the box diagrams shown 
in Fig. 2 in the limit of degenerate mirror leptons. As the box contributions vanish in the 
limit of degenerate mirror quarks, the inclusion of mass splittings in the lepton spectrum 
is a higher order effect. We have numerically verified that it can be neglected for the 
range of mirror fermion masses considered in the analysis. 

Let us begin with the neutral gauge boson contributions. Similarly to AF = 2 transi- 
tions considered in [32,33], only the g^ v part of the gauge boson propagators is relevant. 
The contributions involving k^k u / M^ H cancel each other between the two last sets of 
box diagrams in Fig. 2. Consequently the neutral gauge boson box contributions to X° dd 
and Y° dd are gauge independent. This means that the neutral gauge boson contributions 
to Z^-penguins must vanish in an arbitrary gauge which is confirmed through an explicit 
calculation, as discussed below. 

The result for the box contributions involving turns out to be divergent. As 
box contributions in a renormalizable gauge, like the 't Hooft-Feynman gauge, are finite 
by power counting, the box diagram contributions involving must then be gauge 
dependent. Before giving the result for the unitary gauge calculation, including also 
neutral gauge boson contributions, let us repeat the calculation in the 't Hooft-Feynman 
gauge. 

5.4 Calculation in the 't Hooft-Feynman Gauge 

In the 't Hooft-Feynman gauge also the diagrams with Goldstone bosons have to be 
included. Let us first compute the Z^-penguin diagrams. The O (1) contributions vanish 
as expected and we have to consider 0(v 2 /f 2 ) corrections. Here, as in the unitary 
gauge, there are no contributions from diagrams involving only gauge bosons. On the 
other hand diagrams with Goldstone bosons contribute at 0(v 2 /f 2 ). To this end we 
had to generalize the Feynman rules of [32] to include O (v 2 / f 2 ) corrections to vertices 
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Figure 3: O (v 2 / ' f 2 ) contributions to Zi-penguin in the 't Hooft-Feynman gauge. 

involving Goldstone bosons. It turns out that 0(v 2 /f 2 ) corrections to quark-mirror 
quark-Goldstone boson vertices cancel in the calculation, which implies that the neutral 
gauge boson contributions to the Z^-penguin, not having triple gauge boson vertices 
and corresponding vertices with Goldstone bosons, vanish also in the 't Hooft-Feynman 
gauge. This was to be expected as the box contributions from neutral gauge bosons are 
gauge independent. 

Thus in the 't Hooft-Feynman gauge only two diagrams at O (v 2 /f 2 ), shown in Fig. 3, 
contribute to the Z^-penguin vertex. Using the Feynman rules of Appendix B we find 
that the first diagram in Fig. 3 is divergent with the divergence precisely equal to the 
one found in box diagrams with exchanges calculated in the unitary gauge. 

Including next the finite contributions from the diagrams in Fig. 3 and the finite con- 
tributions from box diagrams with and Goldstone boson exchanges in the 't Hooft- 
Feynman gauge, we confirm the final results for X° dd and y. odd obtained in the unitary 
gauge. 



5.5 Final Results for the T-odd sector 

As described above, we have performed the calculation of the functions X° dd and y. odd 
in the unitary gauge and in the 't Hooft-Feynman gauge obtaining the same results. 
In particular, we found that the left-over divergence obtained in the unitary gauge was 
not an artifact of a non-renormalizable gauge but a physical gauge independent result. 
A similar divergence has been found in the Z^-penguin calculation in the LH model 
without T-parity [28]. We will recall the interpretation of this divergent contribution 
given in [28] in the next section. 

The final results for X° dd and Y° dd in the LHT model are then given as follows: 



X 



odd 



&{r P (z 2 ,y) - r"(z l ,y))+^ , (r p (z 3 ,y) - r P (z u y)) 



y; 



odd 



$\j^(z 2 ,y) - J'*(zi,v))+8>(j'*(z it y) - J^(z u y)) 



(5.1) 
(5.2) 
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where 



J 1 * 



'odd 



ZiS odd + F™( Zi ,y;W H ) 
+a(g( Zi , y; Z H ) + y'\ A H ) + G 2 {z h y; 77)) 

ZiS odd + F^( Zi ,y;W H ) 
a(g{ Zi , y- Z H ) + j/'; A H ) - G 2 (^, y; 77)) 



1 , /i 2 



(5.3) 



Q4P 



(5.4) 



(5.5) 



with the functions F uu , F^, G, G\ and G 2 given in Appendix C and the various variables 
defined as follows 



m 



Hi 



m 



Hi 



M w H 



M Z H ' 



rn 



y 



m 



Hi 



ML' 



a Z j with a 



y =ya, rj 



tan 2 9w 
1 



(5.6) 



(5.7) 



In the unitary gauge the results in (5.1)-(5.4) follow from box diagrams only, since the 
Z^-penguin diagrams do not contribute in this gauge, as discussed in Section 5.2. The 
notation in (5.3) and (5.4) indicates which diagrams contribute to a given function, with 
G2 resulting from diagrams with both Zh and Ah exchanges. In the 't Hooft-Feynman 
gauge the contribution of the Z^-penguin diagram is found to be 

1 v 2 



AJ W = A J"" = AJ 



64 f 



AJ = Zi S odd - 8 Zi R 2 { Z i) + -Zi + 2 Zi F 2 ( Zi ) , 



(5.8) 
(5.9) 



where the functions R 2 and F 2 are given in Appendix C. 

The box diagram contribution involving in the 't Hooft-Feynman gauge can 
simply be obtained from (5.1)-(5.4) and (5.8) using the gauge independence of X° dd and 

\>odd 
1 i 

The formulae (5.1)- (5.4) are the main results of our paper. 



6 The Issue of left-over Singularities 

It may seem surprising that FCNC amplitudes considered in the previous section contain 
residual ultraviolet divergences reflected by the non-cancellation of the 1/e poles at 
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0(v 2 /f 2 ) in our unitary gauge calculation. Indeed due to the GIM mechanism the 
FCNC processes considered here vanish at tree level both in the SM and in the LHT 
model in question. Therefore within the particle content of the low energy representation 
of the LHT model there seems to be no freedom to cancel the left-over divergences as 
the necessary tree level counter terms are absent. 

At first sight one could worry that the remaining divergence is an artifact of the 
unitary gauge calculation. However, an additional calculation in the 't Hooft-Feynman 
gauge convinced us that the found divergence is gauge independent. A similar result has 
been found in the context of the LH model without T-parity in [28] and understood as the 
sensitivity of the decay amplitudes to the UV completion of the LH model. 3 The same 
interpretation can be made here. After all, the LHT model is a non-linear sigma model, 
which is a non-renormalizable description of the low energy behavior of a symmetric 
theory below the scale where the symmetry is dynamically broken. 

We have found explicitly that in the 't Hooft-Feynman gauge the singularity followed 
entirely from the interactions of the Goldstone bosons of the dynamically broken global 
symmetry with the fermions. As the Goldstone bosons in question are the only reminis- 
cences of the spontaneous symmetry breakdown present in the low energy theory, the 
estimate of the size of the divergences through their interactions with fermions in Fig. 3 
should in principle be adequate. However, as emphasized in [28] , the light fermions may 
have a more complex relation to the fundamental fermions of the ultraviolet comple- 
tion of the theory. We refer the reader to [28], where a discussion of this issue and a 
comparison with QCD can be found. 

In what follows we will as in [28] remove 1/e terms from (5.5) and set /i — A to obtain 

1 v 2 A 2 
J Z = JZ = z ^~p lo S mT~ > 

as a minimal estimate of the UV sensitivity of the model. Setting 

A = 4tt/, v = 246 GeV, (6.2) 

we find that for / = 1000 GeV, implying M Wh = 652 GeV, 

JZ = JZ = * ■ 0-006 . (6.3) 

In Fig. 4 we plot J vv and as functions of m Hi for three values of m H i with and 
without J^i v and included. We observe that the divergences constitute a sizable 
fraction of the total result. The coefficient of Zi in the divergent terms J^ v and JZ 

3 A similar singularity has been found independently in [29], in the context of the study of electroweak 
precision constraints. 



21 



J 



. 025 




0.02 



.015 



0.01 



.005 



1000 



m m 



.005 



Figure 4: J vv (lower) and (upper) as functions of mm for values of 

mm =400 GeV(dotted), 500 GeV (solid) and 600 GeV (dashed) with (dark) and with- 
out (light) Jj£ V = J^J included. 

is of the same order of magnitude of the analogous linear coefficient in the convergent 
contributions, but roughly four times larger. Moreover, the linear contribution in the 
range of mirror fermion masses considered is the dominant one, thus explaining the 
important impact of the divergences. At first sight this could imply the loss of the 
predictive power of the theory as our estimate of the divergent contribution is clearly an 
approximation. On the other hand the divergence found has a universal character and 
we can simply write 



and treat 5^ v as a free parameter. Assuming that 5di v encloses all effects coming from the 
UV completion, which is true if light fermions do not have a more complex relation to the 
fundamental fermions of the UV completion that could spoil its flavour independence, 
one can in principle fit S^iv to the data and trade it for one observable. At present this 
is not feasible, but could become realistic when more data for FCNC processes will be 
available. 

On the other hand, implementing T-parity removes all divergences from the T-even 
sector. This is easy to understand. The only new T-even particle is T + which can be 
thought of as an arbitrary singlet field mixing with the SM top quark, independently 
of the non-linear sigma model. Of the "pion" matrix EE, only the SM Higgs doublet is 
present in the T-even sector, and all modifications in its couplings appear due to the 
mixing of T + with t. Thus the T-even sector of the LHT model is effectively decoupled 
from the breaking 577(5) — > SO (5) of the non-linear sigma model, which has been the 
basic reason for the appearance of the singularity described above and in [28]. 



J! 



ruu _ tmm 
div ~~ u div 




(6.4) 
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7 B X s £ + £- 



7.1 Preliminaries 

The branching ratio for the rare decay B — > X s £ + £~ depends in the SM on the functions 
Ysm, Zsm and D' su with the latter relevant also for the B — > X s 7 decay. The formulae 
for the branching ratio are very complicated and will not be presented here. They can be 
found in [57], where also the formulae for the forward-backward asymmetries are given. 
In the LHT model the function Ysm is generalized to Y s calculated in the previous section, 
whereas -Dlht ^ as Deen calculated in [33] and is given for completeness in Appendix C. 

What remains to be calculated is the function Z s that we defined in (3.5). The SM 
contribution can be written as 

Z SM = C + ^D (7.1) 

with C and D given in the 't Hooft-Feynman gauge in Appendix C. Z SM is gauge 
independent. We will also need the QCD-penguin function E that can be found in 
Appendix C as well. 



7.2 T-even Sector 

In the case of the T-even sector it is useful to work in the unitary gauge. The function 
Ceven can be extracted from [28] by imposing T-parity with the result 

9 9/ 9 
X L V n ( X t — Xt V 



/-(even _ _L ^_ q Zl A _____ 

^unitary — g p ° cvcn I ^ U 2^T j 2 

x\ v 2 f — 6 — 5x t + 5x 2 — 3xt + 3x t xr 



where 



16 P V 2(x t - 1) 

8x t — lOxf + 5rr? 
H ^ _ -g 2 log x t - [Ax t + x T ) log x T 

x 2 v 4 / 3 \ 
+ YJI XT ( ~2 d2 + x l + d 2 lQ g^ J ' ( 7 - 2 ) 



1 /x 2 

^cvcn = - + log , (7.3) 



and d 2 has been defined in (2.34). 

The function -D e ven in the LHT model can be obtained with the help of D SM in the 
unitary gauge. To our knowledge the latter function has never been given in this gauge 
in the literature. It can be found with the help of [28] as follows. 
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From the gauge independence of Z$m we know that 



Csm — Co + -Qsm , D SM — D — 2g SM , 



(7.4) 



where qsm is gauge dependent and vanishes in the 't Hooft-Feynman gauge. It has been 
calculated in the unitary gauge in [28] with the result 



1 -3x 2 + 17x t 8x 2 - x 3 

8*^ ven 



16(1 -x t ) 8(^-l) 2 
Consequently, using _D in Appendix C, we find 

-153xt + 383^ _ 2A5x 3 + 27 x f 



(7.5) 



nSM / \ _ ^io 

^unitary K^t J — ^ ^evcn T 



72(x t - l) 3 



16 - 64a; t + 36x 2 t + 93x 3 t - 8Axf + 9x 5 t 
36{x t - l) 4 



logs*. 



(7.6) 



Proceeding then as in the case of B — > X s ^ in [28] but including also f 4 // 4 corrections 
to the diagrams with internal T + exchanges we find 



D 



unitary 



/ 



2 X L 



1 +2d 2 j^ ] D 'unwary (x T ) - -Dunftary ( X t) 



(7.7) 



and, dropping 0{v / f ) terms, 



7-)even ^ 2 

^unitary f 2 



^ (l + 2<i 2 ^ ) S even - Z^L, y 0. 



+^1 



41 - 24 log x T , x T 
~8 
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f ( l + 2rf 2 — ) (3-21ogx T ) 



and subsequently 



^evcn — "^unitary ^-^unitary 



(7.8) 



(7.9) 



which is gauge independent. The divergence S even cancels in (7.9) so that Z even is finite, 
in agreement with our statement in Section 6. 

7.3 T-odd Sector 

In the T-odd sector it is useful to work in the 't Hooft-Feynman gauge. Let us denote 



Zodd(zi) — Codd(^j) + ^-Dodd(^i) , 
then, from (5.8) and (5.9), we find in the 't Hooft-Feynman gauge 



c M = aj- = 1^ 



ZiS odd - 8ziR 2 {zi) + -Zi + 2ziF 2 {zi 



(7.10) 



(7.11) 
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with S dd defined in (5.5) and the functions R2 and F2 given in Appendix C. 

D dd{ z i) is then found in analogy to our calculation of the B — > X s ^ decay in 
We obtain 

D dd{Zi 



D (zi) - ^E ( Zl ) - ^o(^) 



with D and E given in Appendix C. Finally we have 



7odd 



£.2 ^ (^odd(^2) — Z od d(Zl)) + £3^ (Z dd{Zz) — ^odd(^l)) 



(7.12) 



(7.13) 



As D Q dd(zi) is finite, the divergence in C Q dd(zi) remains in Z odd . Its estimate can be done 
as outlined in Section 6. 



8 Kl — > tt°£ + £- 

The rare decays Kl — > n e + e~ and Kl — > 7r°/x + / u _ are dominated by CP-violating 
contributions. In the SM the main contribution comes from the indirect (mixing-induced) 
CP violation and its interference with the direct CP- violating contribution [58-61]. The 
direct CP-violating contribution to the branching ratio is in the ballpark of 4 ■ 10~ 12 , 
while the CP conserving contribution is at most 3 • 1CT 12 . Among the rare K meson 
decays, the decays in question belong to the theoretically cleanest, but certainly cannot 
compete with the K — > txvv decays. Moreover, the dominant indirect CP-violating 
contributions are practically determined by the measured decays K$ — > tt°£ + £~ and the 
parameter ek- Consequently they are not as sensitive as the Kl — > ts^vv decay to new 
physics contributions, present only in the subleading direct CP violation. However, as 
pointed out in [37], in the presence of large new CP- violating phases the direct CP- 
violating contribution can become the dominant contribution and the branching ratios 
for K L — > ir°l + l~ can be enhanced by a factor of 2-3, with a stronger effect in the case 
of Kl ->■ ttW" [60,61]. 

Adapting the formulae in [59-62] with the help of [37] to the LHT model we find 

Br (K L n°l + r) = {C e dil ± C( nt \a s \ + Cf nix |a s | 2 + C e CPC ) • 10" 12 , (8.1) 

where 

C e dh = (4.62 ± 0.24)(w 7 V + oo 2 7A ) , = (1.09 ± 0.05) (w^ + 2.32^) , (8.2) 

Cf nt = (11.3 ± 0.3)u 7V , Cf nt = (2.63 ± 0.0Q)u 7V , (8.3) 

C e mix = 14.5 ± 0.05 , C^ lx = 3.36 ± 0.20 , (8.4) 

C£ PC ~0, Cg PC = 5.2±1.6, (8.5) 

|a s | = 1.2 ± 0.2 (8.6) 
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with 



U7V ~ 27 



bJ-jA = 



P +~ 



\Y, 



it- 



sin (3- 



sin 2 9 W sin(/3 - (3 S 



A\Z 



sin (3 



K 



sin(/3 - (3 S ) 



Im A t 



1.4 • 10" 4 



\Y, 



K 



sin /3 



K 



2ir sin 2 6>vf sin(/3 - f3 s 



Im A t 



1.4- 10" 4 



where P = 2.88 ± 0.06 [63] includes NLO QCD corrections and 



K 



9 Z 



(8.7) 
(8.8) 

(8.9) 



tf = p-p.-e$, P§ = p-P. 

with Z K defined in (3.5) and obtained from Z s by changing ^ to 

The effect of the new physics contributions is mainly felt in uja, as the corresponding 

contributions in uj-jy cancel each other to a large extent. 
The present experimental bounds 



Br[K L ->■ 7T°e + e~) < 28 • lO" 11 [64], Br(K L ->■ ttVV") < 38 • 10" 

are still by one order of magnitude larger than the SM predictions [62] 
Bt(Kl 7r e + e -)sM = 3.54±g:H (l-56±g : 2) • lO" 11 , 
Bt(Kl ttV/Osm = 1.4llg;| (0.95^; 22 ) • lO" 11 

with the values in parentheses corresponding to the "— " sign in (8.1). 



ii 



[65] 
(8.10) 

(8.11) 
(8.12) 



9 Benchmark Scenarios for New Parameters 
9.1 Preliminaries 

In what follows, we will consider as in [33] several scenarios for the structure of the Vna 
matrix and the mass spectrum of mirror fermions with the hope to gain a global view 
about the possible signatures of mirror fermions in the processes considered and of T + 
present in the T-even contributions. 

In the most interesting scenario considered in [33] (Scenario 4 below), the mixing 
matrix Vud differed significantly from Vckm- It could have a large non- vanishing complex 
phase 5f 3 , while the phases Sf 2 and ^3, with smaller phenomeno logical impact, were set 
to zero. In this scenario large CP-violating effects in B s decays have been found. In 
particular, the CP asymmetries and A$ L could be enhanced by an order of magnitude 
with respect to the SM expectations. 

In the next section we will be primarily interested in calculating the observables 
considered in the previous sections in the scenarios defined in [33] . In particular, it will be 
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interesting to see how the CMFV correlations between K°, B® and B® systems [40,67] are 
modified when new sources of flavour and CP violation are present. The parameterization 
of various decays in terms of the functions Xi, Yi and Zi that we defined in Section 3 is 
very useful for such tests. 

The main purpose of our numerical analysis is to have a closer look at six scenarios, 
with the first five considered already in our previous study of particle-antiparticle mixing 
and B — y X s ^. We will recall these five scenarios and introduce a sixth one which has 
particularly interesting effects in K — y nvv and K L —y ir°£ + l~ decays. In all these six 
scenarios we set the phases Sf 2 and 8 23 to zero. The values of the observables that can 
only be produced by allowing these phases to vary freely, are covered by our general scan 
anyway. This simplification therefore does not restrict the generality of the analysis. We 
will see that Scenarios 4 and 6 turn out to be the most appealing, since they respectively 
provide large enhancements in the B s and K systems. Spectacular effects in both B and 
K systems, however, are not simultaneously allowed in a single scenario. Therefore, we 
complete the numerical analysis with a general scan over mirror fermion masses and Vna 
parameters, with also Sf 2 and S 23 free to differ from zero, finding some interesting points 
where significant enhancements in both B and K systems occur. 

9.2 Different Scenarios 

Here we just list the scenarios in question: 
Scenario 1: 

In this scenario the mirror fermions will be degenerate in mass 

m m = m H2 = m m (9.1) 

and only the T-even sector will contribute. This is the MFV limit of the LHT model. 
Scenario 2: 

In this scenario the mirror fermions are not degenerate in mass and 

V Hd = Vckm ■ (9.2) 

In this case there are no contributions of mirror fermions to D° — D° mixing and flavour 
violating D meson decays, and 

g> = A ( 9 ) f = A («) ) (9 . 3) 

with q = d,s and no index q in the K system. As discussed in [33] and below this 
scenario differs from the MFV case. 
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Scenario 3: 

In this scenario we will choose a linear spectrum for mirror quarks 



m m = 400 GeV, m H2 = 500 GeV, m H3 = 600 GeV , (9.4) 

set m m = 500 GeV and take an arbitrary matrix Vna but with the phases Sf 2 and 5 23 
set to zero. We stress that similar results are obtained by changing the values above by 
±30 GeV, with similar comments applying to (9.5) below. In the remaining scenarios we 
will modify the mirror quark spectrum but keep m m = 500 GeV. 

Scenario 4: 

This was our favorite scenario in which large departures from the SM and MFV in B s 
decays could be obtained and some problems addressed in [33] could be solved, with small 
effects in the experimentally well measured quantities AMk and Sk- In this scenario 

m m « m H2 = 500 GeV , m H3 = 1000 GeV , (9.5) 

-^<s? 2 <0.99, 5- 10~ 5 < s^ 3 < 2- 10" 4 , 4- 10~ 2 < s( 3 < 0.6, (9.6) 
v2 

Sf 2 and 5 23 are set to zero, while the phase Sf 3 is arbitrary. The hierarchical structure of 
the CKM matrix 

sia < s 23 < sia , (CKM) (9.7) 

is changed in this scenario to 

4«4<4, (v Hd ) (9.8) 

so that Vhcl looks as follows: 

/ c 12 b \2 A 13 e \ 



V U d = 



s \2 u 12 ' 5 23 

\-cf 2 s d 13 e i5 ri -sf 2 sj 3 e i5 ™ 1 / 



(9.9) 



We would like to stress that with the degeneracy m H1 m^ 2 the T-odd contributions in 
ek proportional to Im(£ 2 ) and Re(^) vanish, and only the T-odd term proportional to 
Im(^ 3 )Re(^ 3 ) contributes. Being Im(£ 3 ) = sf 3 c 23 s 23 sm8f 3 , the hierarchy chosen in this 
scenario for Vnd, with s 23 1, has the advantage of suppressing mirror fermion effects 
in e K , allowing at the same time large CP- violating effects in the — B° s system [33]. 
Furthermore AM S can be smaller than its SM value in this scenario, and interesting 
effects in the B® — B° d system are also found. 

It will be interesting to see whether in this scenario large departures from the SM 
expectations for rare decays can be obtained. 
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Scenario 5: 

In all the previous scenarios we will choose the first solution for the angle 7 from tree 
level decays as given in (10.1) below so that only small departures from the SM in the 
B% — B^ system will be consistent with the data. In the present scenario one assumes the 
second solution for 7 in (10.1) in contradiction with the SM and MFV. We have shown 
in [33] that the presence of new flavour violating interactions could still bring the theory 
to roughly agree with the available data, in particular with the asymmetry S^k s - In 
spite of that, the combined measurements on A$ L and cos(2/3 + 2(pB d ) and the indirect 
experimental estimate of A$ L make this scenario very unlikely [33,68], such that we will 
not consider it any further. 

Scenario 6: 

In studying this scenario we aim to enhance mirror fermion contributions to rare K 
decays, keeping negligible effects in the experimentally well measured quantities AMk 
and Sk- To this purpose we choose the mirror fermion masses as in Scenario 4 (see (9.5)) 
since the near degeneracy between mni and m,H2 helps to suppress mirror fermion effects 
in AM K . 

Concerning e k , we recall that with the degeneracy mm ~ m H2 the T-odd contribu- 
tions proportional to Im(£ 2 ) and Re(^) vanish, and only the T-odd term proportional 
to Im(£ 3 )Re(£3) contributes. In Scenario 4 the hierarchical structure of Vna is chosen as 
to satisfy Im(£ 3 ) ~ 0. Here in Scenario 6, instead, we suppress mirror fermion effects in 
Ek due to the second and third generations, by requiring Re(£ 3 ) = 0. Setting also in this 
scenario the phases Sf 2 and 5 23 to zero, the explicit expression of the real part reads 



Re(&) = -44 (4 2 - 4 3 2 s d 13 2 ) + (cos s d 13 ) 444 (4 2 - 4 2 ) , (9.10) 



We note that while the value of 6f 2 is fixed to 45° by (9.11), 6 23 and 8f 3 have no 
specified value nor order of magnitude, but (9.12) implies that only one of them is a free 
parameter. The matrix Vud can then be expressed in terms of the two free parameters 



which vanishes for 8f 2 , 6 23 and 9f 3 (chosen in the first quadrant) satisfying 




(9.11) 



d _ °13 



(9.12) 
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9f 3 and 8f 3 as 



/ 



V, 



Hd 



V2 



;i + 4 2 e*) 



V2 

(1 " sf 3 2 



s" 3 e 13 



c 13 



(9.13) 



/ 



Its structure becomes much simpler if the angle 9f 3 is sufficiently small, i.e., sf 3 < 0.1, 



and reads 



V, 



Hd 



1 
V2 



1 
1 



s 13 e 



A 



'13 



\V2 



g*^13 



) 



(9.14) 



-3|(l + e *&) 1 y 

As we will see in Section 10 the very different structure of Vud when compared with 
Vckm implies enhancements in rare K decays, without introducing problematic effects 
in AM K and e K . Moreover, as V Hd in (9.14) has a different structure also from the (9.9) 
one of Scenario 4, the new physics effects in the — and mainly in the B° s — B° s 
system, turn out to be small although visible. 

General Scan: 

As shown in Section 10, Scenarios 4 and 6 turn out to be the most interesting ones 
with, respectively, large new physics effects in the B s and K systems. Such visible 
enhancements follow from the structure of Vnd, primarily required to satisfy the Ek and 
AM K constraints, through Im(£ 3 ) ps in Scenario 4 and through Re(£ 3 ) = in Scenario 
6. A further consequence of the Vna structure is that in Scenario 4 spectacular effects 
can be obtained in the B s system but not in the K system and vice versa in Scenario 
6. An even more interesting picture would be the simultaneous manifestation of large 
enhancements in both B and K observables. In order not to miss such a possibility, in 
addition to the scenarios described above, we have performed a general scan over mirror 
fermion masses and V Hd parameters. To have a global view of the most general LHT 
effects, we have allowed here the phases Sf 2 and 8% 3 to differ from zero. Qualitatively 
their effect is not significant, although they can help in achieving very large effects in 
certain observables. We find that there exist some sets of masses and V Hd parameters 
where the new physics effects turn out to be spectacular in both B and K systems. 
We note that they do not really constitute a scenario, they rather appear in the plots 
shown in the next section as isolated (blue) points. In contrast to previous scenarios, in 
fact, the blue points corresponding to large new physics effects are quite sensitive to the 
particular configuration of mirror fermion masses and Vnd parameters. 
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10 Numerical Analysis 



10.1 Preliminaries 

In our numerical analysis we will set | V us | , | V c b | and | V u b | to their central values measured 
in tree level decays [69,70] and collected in Table 1. 



\v ub 


= 0.00429(29) 




G F = 1.16637(1) • 10" 5 


GeV" 2 


\Vcb\ 


= 0.0416(9) 


[69] 


M w = 80.425(38) GeV 




A = 


\V US \ = 0.225(1) [70] 


a = 1/127.9 




\Vts\ 


= 0.0409(9) 


[68] 


sin 2 6 W = 0.23120(15) 




A = 


0.822(16) 




m K o = 497.65(2) MeV 




Rb = 


= 0.447(31) 




m Bd = 5.2794(5) GeV 




P = 


26.3(21)° 




m Bs = 5.370(2) GeV 




Ps = 


-1.28(7)° 




F K = 160(1) MeV 


[71] 


m c = 


= 1.30(5) GeV 




F Bd = 189(27) MeV 




m t = 


= 163.8(32) GeV 




F Bs = 230(30) MeV 


[72] 



Table 1: Values of the experimental and theoretical quantities used as input parameters. 

As the fourth parameter we will choose the angle 7 of the standard unitarity triangle 
that to an excellent approximation equals the phase #ckm in the CKM matrix. The angle 
7 has been extracted from B — > D^K decays without the influence of new physics with 
the result [68] 

7 = (71 ±16)°, 7 = -(109 ± 16)°. (10.1) 

Only the first solution agrees with the SM analysis of the unitarity triangle, while the 
consistency of the second solution with data has been investigated within Scenario 5 in 
our previous LHT analysis [33]. It turns out that the combined measurements on A§ L 
and cos(2/3 + 2ip Bd ) and the indirect experimental estimate of A$ L make this scenario 
very unlikely [33,68]. 

We will consider here only the first solution, whose uncertainty is sufficiently large 
to allow for significant contributions from new physics. The value of f3 in Table 1 is 
obtained from and the first solution of 7, i.e. from tree level decays only and is not 
affected by an eventual new physics phase. Its difference from the value of (3 obtained 
from the S^k s asymmetry, (5(^Ks) = 21.2 ±1.0, constitutes the "sin 2/3 problem" which 
can be solved only in Scenarios 3 — 6 [33] . 

For the non-perturbative parameters entering the analysis of particle-antiparticle 
mixing we choose and collect in Table 1 their lattice averages given in [72], which combine 



31 



unquenched results obtained with different lattice actions. Other parameters relevant 
for particle-antiparticle mixing can be found in [33]. 

In order to simplify our numerical analysis we will, as in [33], set all non-perturbative 
parameters to their central values and instead we will allow AMk, £k, AM d , AM s and 
S i>Ks to differ from their experimental values by ±50%, ±40%, ±40%, ±40% and ±8%, 
respectively. In the case of AM s /AM d we will choose ±20% as the error on the relevant 
parameter, £, is smaller than in the case of AM d and AM S separately. The relevant 
expressions are given in [33]. These uncertainties could appear rather conservative, but 
we do not want to miss any interesting effect by choosing too optimistic non-perturbative 
uncertainties. The constraints from B — > X s ^ and B — > X s l + £~ are also taken into 
account. They turn out to be easily satisfied, within the present uncertainties, and 
therefore to have only a minor impact. 

In Scenarios 3 — 6, the parameters / and x L will be fixed to / = 1000 GeV and 
xl = 0.5 in accordance with electroweak precision tests [30]. Varying the breaking scale 
/ would obviously modify our results. The effect, however, turns out to be much smaller 
than one would naively expect from the v 2 If 2 — "scaling". In other words, lower values 
of / do not allow arbitrarily large NP contributions, since the constraints imposed from 
the available data become more stringent in this case. 

10.2 The MFV Scenario 1 

Let us consider first the case of totally degenerate mirror fermions. In this case the odd 
contributions vanish due to the GIM mechanism [47], the only new particle contributing 
is T + and the LHT model in this limit belongs to the class of MFV models. As only 
the T-even sector contributes, the new contributions to all FCNC processes are entirely 
dependent on only two parameters 

xl, f- (10.2) 
Moreover, all the dependence on new physics contributions is encoded in the functions 

Veven = -<^SM ± -Y even , ^cven = ^SM ± ^even j Z even = 

(10.3) 

There exist strong correlations between various processes that are characteristic for 
models with MFV. 

It should be emphasized that in this scenario the "sin 2/3 problem" cannot be solved 
as it is a MFV scenario and that AM S > (AM S ) SM , which is not favored by the CDF 
measurement [73], as well as AM d > (AM d ) SM . In [74] the relations AM s4 > (AM s4 ) SM 
have been proven to be valid in constrained MFV, where flavour violation is governed 
entirely by the Yukawa interactions and there are no new operators beyond the SM ones, 
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Figure 5: X even /X SM (left) and Y even /Y SM (right) as functions of xl, for various values 
of f = 1, 1.2, 1.5 and 2TeV from top to bottom. The bands underlying the curves show 
the allowed ranges after applying electroweak precision constraints [30]. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

x L x L 



Figure 6: Br(K L n°uu)/Br{K L n vv) SM (left) and Br(K + n + uu)/Br(K + -> 
7T + b>i?)sM (right) in Scenario 1, as functions of xl for different values of f = 1,1.2,1.5 
and 2 TeV from top to bottom. The bands underlying the curves show the allowed ranges 
after applying electroweak precision constraints [30]. 

and, therefore, have been expected for this scenario. We specify that in the numerical 
analysis of this scenario the S^k s constraint is left out while the AM^ S ones are taken 
into account. 

In Fig. 5 we plot X eYen / X^m. and F even /YsM as functions of xl for various values of 
/. We observe that the new physics contributions amount to modifications of the SM 
functions by at most 9% and 14%, respectively, and of the corresponding rare decay 
branching ratios by at most 18% and 28%. As an example we show in Fig. 6 Br{Ki — > 
ti°i/V) / Br(Ki — > tt uu)sm an d Br(K + — > -k + vv)/ Br(K + — > 7t + z/z/)sm as functions of xl 
for different values of /. We observe that slightly larger effects can occur in BriK^ — > 
^vv) relative to Br(K + — > tt + vv). 
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10.3 Scenario 2 



The new contributions to FCNC processes are in this scenario entirely dependent on 
only six parameters 

xl, f, rn H i, m H 2, m H3 , m m , (10.4) 

in addition to m t and the CKM parameters that we set to the central values obtained 
from tree level decays. 

In spite of the fact that in this scenario V Hd = Vqkm, it does not belong to the class 
of MFV models. The point is that breaking the degeneracy of mirror fermion masses 
introduces a new source of flavour violation that has nothing to do with the top Yukawa 
couplings. Only if accidentally the contributions proportional to ^ = X[ 9 ^ dominate 
the new physics contributions, one would again end up with a scenario that effectively 
looks like MFV. However, as the mirror spectrum can be generally different from the 
quark spectrum and not as hierarchical as the latter one, the terms involving ffi = ^ 
in the formulae of the previous sections cannot be neglected, while this can be done in 
the T-even contributions. Moreover, as Ac 9 " 1 are different from X[ q \ that dominate the 
SM contributions, even in this simple scenario the usual MFV relations between K, B d 
and B s systems will be violated. 

Specifically, for q = d, s, X c are of the same order of mag nitude as X { t q) and the MFV 
relations between B d and B s systems turn out to be only weakly violated. On the other 
hand, X^/\[ K ^ = 0(4 ■ 10 2 ), implying that not only the MFV relations between K and 
Bd, s systems are strongly violated, but also the rate Br(K + — > -k + vv) can be significantly 
enhanced in this scenario, relative to the SM and Scenario 1. In Br(K + — > 7i + uu), in 
fact, the T-odd contribution proportional to Re(Ai^) can have a significant effect, since 
Re(Ai^) is much larger than lm(X[ K ^) and Re(A^). In BriK^ — > tx®vv\ instead, only 
the imaginary part of the CKM contributions enters and, since 
the T-odd contribution can only yield a slight suppression. 

As an example, we show in Fig. 7 Br(K L — > -K vis)/Br(K L — > 7t z/z/) S m and Br(K + — > 
tt + uP)/ 'Br(K + — > 7r + z/z/) SM as functions of x L , choosing / = 1 TeV and scanning over 
mirror fermion masses. We point out that the central values of the SM predictions 
appearing in the ratios shown in Fig. 7 differ from those quoted in (3.13) and read 
Br{K + -> vr+z/z/) = 8.7 • 1CT 11 , Br{K L n°uu) = 4.1 • HT 11 . This difference comes 
from the CKM inputs that in the present analysis are taken from tree-level decays only. 
We observe that Br(K L — > rPvv) can be enhanced at most by 17% relative to the SM 
prediction, with stronger new physics effects at higher values of the Xl parameter. In 
Br(K + — > it + vu), instead, no clear dependence on the xl parameter can be seen, while 
larger (of a factor 5) enhancements as well as suppressions of an order of magnitude can 
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Figure 7: Br(K L 7r°vv)/Br(K L n°vv) SM (left) and Br{K + 7i + uu)/Br(K+ -> 
7r + z/z/)5M (right) in Scenario 2, as functions of xl, choosing f = 1 TeV and scanning over 
mirror fermion masses. The solid line in the right plot represents the SM prediction and 
separates the regions where Br(K + — > 7t + vu) is suppressed or enhanced relative to the 
SM. 

be obtained. 

10.4 Breakdown of the Universality 

In MFV models the functions X{, Yj and Zj are independent of the index i — K,d, s. 
Consequently, they are universal quantities implying strong correlations between ob- 
servables in K, Bd and B s systems. The presence of mirror fermions in the LHT model 
generally breaks this universality, as we have already seen in Scenario 2. 

In Fig. 8 we show the ranges allowed in different scenarios in the space (|X 8 |, 
Here and in all the following plots, Scenarios 3, 4 and 6 are represented by red, green 
and brown points, respectively, while blue points stand for the general scan. The solid 
line represents the MFV relation \X S \ — \Xk\, with the black point corresponding to the 
SM prediction Xsm = 1-49 and the light blue point showing, for illustrative purposes, 
the Scenario 1 result. The departure from the solid line gives the size of non-MFV 
contributions allowed in the various, differently coloured, scenarios. We observe that 
roughly 

1.40 < \X S \ < 1.75, 0.7 < \X K \ < 4.7, (10.5) 

implying that the CP-conserving effects in the K system can be much larger than in 
the B s system. A very similar behavior is found for the functions Yi and Zj with larger 
effects in the K system relative to the B^ s systems. For instance 

0.89 < m| < 1.17, 0.41 < \Y K \ < 3.9 , (10.6) 

to be compared with |y s | = \Yk\ = 0.95 in the SM. 
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In Fig. 9, then, we show the allowed ranges in the space which roughly 

turn out to be 

-3.5° < 0^ < 3.5° , -130° < 6£ < 55° , (10.7) 

implying that the CP-violating effects in the b — > s transitions are tiny, while those 
in Kl decays can be very large. An analogous pattern is found for the Y Sy x and Z St K 
functions. For the absolute values, the ranges in the Bd and B s systems are similar, 
while in the Bd system larger values for the phase (±13°) can be reached. 

From the discussion of the previous plots it is evident that mirror fermion contribu- 
tions break universality, and in a scenario-dependent way. Furthermore, we would like 
to stress and explain the origin of larger effects found in the K system relative to the 
Bd, s systems. Looking at the X« expression in (3.3) one sees that the mirror fermion 
contribution is enhanced by a factor As X[ K ^ ~ 4 • 10~ 4 , whereas A^ ~ 1 • 1CT 2 

and A^ ~ 4 ■ 1CT 2 , we naively expect the deviation from X SM in the K system to be by 
more than an order of magnitude larger than in the B d system, and even by two orders 
of magnitude larger than in the B s system. Analogous statements are valid for the Yi 
and Zi functions. 

In view of the smallness of the new physics contributions in b — > s transitions it is 
easy to satisfy the constraints from B — > X s ^ and B — > X s £ + £~ that turn out to be close 
to the SM expectations. Therefore we will not further discuss them. 

10.5 The K <kvv System 

In Fig. 10 we show the correlation between Br(K + — > 7i + uu) and BriK^ — > tx^vv) for 
the scenarios in question. The experimental la-range for Br(K + — > 7i + uu) [52] and 
the model-independent Grossman-Nir (GN) bound [75] are also shown. We observe that 
even for the most general case, there are two branches of possible points. The first one is 
parallel to the GN-bound and leads to possible huge enhancements in Br(K^ — > 7r° vv) 
so that values as high as 5 • 10~ 10 are possible, being at the same time consistent with 
the measured value for Br(K + — > 7i + uu). Within Scenario 6 (brown points), this branch 
reduces to a straight line. The second branch corresponds to values for Br(K L — > ix^vv) 
being rather close to its SM prediction, while Br(K + — > tt^uu) is allowed to vary in the 
range [1 • 10~ n , 5 • 10~ 10 ], however, values above 4 • 10~ 10 are experimentally not favored. 
We note that within Scenario 4 (green points), Br(K L — > n^vv) is fixed to the T-even 
contribution and close to the SM value, and Br(K + — > is always smaller than in 

the SM so that the GN-bound can be reached. 

In Fig. 11 we show the ratio Br(K L — > TC vis)/Br(K + — > -n+vv) as a function of the 
phase $xi displaying again the GN-bound. We observe that the ratio can be significantly 
different from the SM prediction, with a possible enhancement of an order of magnitude. 
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Figure 10: Bt{Kl — > 7r°z/z/) as a function of Br(K + — > tt + vv). The shaded area repre- 
sents the experimental la -range for Br (K + — > 7i + uu). The GN-bound is displayed by the 
dotted line, while the solid line separates the two areas where BriK^ — > 7i uu) is larger 
or smaller than Br(K + — > 7i + uu). 



Br (K L 71° v v) Br (K L ^ 7T° v v) 

Br (K + 7T + v v) Br (K+ ^ 7T+ V v) 




Figure 11: Br(Ki —> ti°vv) / Br{K + — > ti + vv) as a function of fix > ^ n general scan 
and Scenarios 3 ; 4 and 6, respectively. The dashed line represents the GN-bound. 
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The two branches of Fig. 10 can be distinguished also here. In particular, points gener- 
ated in Scenario 6 appear only in the left branch and can not reach the GN-bound, while 
points belonging to Scenario 4 populate the right branch and approach this bound. 
The most interesting implications of this analysis are: 

• If Br(K + — > 7T + z/z/) is found sufficiently above the SM prediction but below 2.3 • 
10~ 10 , basically only two values for BriK^ — > 7T uu) are possible within the LHT 
model. One of these values is very close to the SM value in (3.13) and the second 
much larger. 

• If Br(K + — > is found above 2.3 • 1CT 10 , then only Br(K L — > ifivv) with a 
value close to the SM one in (3.13) is possible. 

• If Br(K + — > Ti + uu) is found above the SM value, Scenario 4 will be ruled out. 

As Scenario 4 was our favorite scenario in the analysis of AF = 2 processes in [33] , with 
spectacular new physics effects in and A$ L asymmetries, let us next have a closer 
look at the correlations between and the K — > ixvv decays. 

10.6 S^^ and K — > ttuu 

In Figs. 12 and 13 we show the correlation between and BriK^ — > 7i uu) and 
Br(K + — > 7i + uu), respectively. We observe that in Scenario 4 (green points), in which 
S,^^ can be significantly enhanced, BriK^ — > 7T uu) is very close to the SM value, while 
Br(K + — > 7i + uu) is suppressed as already found previously. On the other hand in 
Scenario 6 (brown points), both Br(K^ — > 7r°z/z/) and Br(K + — > 7i + uu) can be strongly 
enhanced, while is very close to the SM value. Only the general scan (blue points) 
can yield a simultaneous enhancement of these three observables. In order to see the 
triple correlation in question even better for all the scenarios, we show in Fig. 14 only 
those points of Fig. 10 for which > 0.1. It is evident that it is rather difficult to 
obtain simultaneously large values of the three observables in question. Still some sets 
of parameters belonging to the general scan exist for which this is possible. 

10.7 B ->> X^ d vi? 

From the discussion of Section 10.4 we conclude that the branching ratios for the B — > 
X S)d vv decays can be enhanced by at most 35% over the SM predictions. Moreover, we 
find that 

2 

< 1.56 , (10.8) 
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Figure 12: Bt^Kl — > 7r°z/z/) as a function of S^, in the general scan and Scenarios 3, 4 
and 6, respectively. 
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Figure 13: Br(K + — > ti + vv) as a function of S^. The shaded area represents the 
experimental la-range for Br(K + — > n + uu). 
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Figure 14: Br{Ki — > TT ui?) as a function of Br(K + — > 7i + uu), showing only the points 
that satisfy > 0.1. Like in Fig. 10, the shaded area represents the experimental 
la-range for Br(K + — > tt+vv), the GN-bound is displayed by the dotted line, while 
the solid line separates the two areas where BriK^ — > n^vv) is larger or smaller than 
Br(K + tt+pu). 

implying that the MFV relation between the ratio of the branching ratios in question 
and the CKM parameters given in (3.25) can be significantly violated. 

10.8 Bd,s versus K — > Ttvi> 

We next investigate possible correlations between Bd. s — > ari d K ~^ kvv. In 

particular, we show in Fig. 15 the first correlation that will be accessible to future 
experiments: Br(B s — > fi + fi~)/Br(B s — > )sm as a function of Br(K + — > 7i + uu). 
The experimental la-range for Br(K + — > tt + uu) [52] is represented by the shaded area 
and the SM prediction by the dark point. Br(B s — > can be about 50% larger than 

in the SM, while more pronounced effects are possible in Br(K + — > 7i + uu). Scenarios 
4 (green points) and 6 (brown points) turn out to be again distinguishable through 
this correlation. As expected from our previous discussion, Scenario 4 allows larger 
effects in the B s system, while Scenario 6 is characterized by significant (of a factor 5) 
enhancements in the K system. 
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Br (B a \T ) 




Figure 15: Br(B s — > / u + ^r) / 'Br(B s — > /j> + ^~)sm as a function of Br(K + — > 7T + uu). 
The shaded area represents the experimental la-range for Br(K + — > 7i + uu) and the dark 
point shows the SM prediction. 




Figure 16: Br(Ki — > 7T°fi + {i ) as a function of BriK^ — > 7r°e + e ). 
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Figure 17: Br(Ki — > 7r°e + e _ ) (upper curve) and Br(K^ — > 7r°fi + fi~) (lower curve) as 
functions of Br(K L — > n^vv). The corresponding SM predictions are represented by dark 
points. 

10.9 The K L ir°£ + £- System 

In Fig. 16 we show the correlation between BriK^ — > ir°e + e~) and BriK^ — > 7r°yU + /i~) 
that has been first investigated in [60,61] in the framework of [37]. This correlation is 
only moderately sensitive to the scenario considered. We observe that both branching 
ratios can be enhanced up to a factor two, over the SM values (black point) in (8.11) 
and (8.12). 

10.10 K, ir°£ + £- versus K L ttW 

In Fig. 17 we show Br(Ki — > 7r°e + e - ) and BriK^ — > 7r°/i + / u _ ) versus BriK^ — > tt vv). 
We observe a strong correlation between — > 7i°£ + £~ and Ki — > ifivv decays that we 
expect to be valid beyond the LHT model, at least in models with the same operators 
present as in the SM. We note that a large enhancement of BriK^ — > 7r°z/p) automat- 
ically implies significant enhancements of Br(Ki — > 7r°£ + £~) and that different models 
and their parameter sets can than be distinguished by the position on the correlation 
curve. Moreover, measuring BriK^ — > tt°£ + £~) should allow a rather precise prediction 
of Br(KL 7r°z/p) at least in models with the same operators as the SM. 
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Figure 18: The ratio r as a function of 5 
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Figure 19: sin 2/3^/ sin(2/3 + 2(p Bd ) as a function of 5f 3 . 
10.11 Violation of Golden MFV Relations 

In Fig. 18 we show the ratio r of (3.20) as a function of 5f 3 . Its departure from unity 
measures the violation of the golden MFV relation between Bd )S — > H + fJ>~ decays and 
AM^ S in (3.20). We observe that r can vary in the range 



0.6 < r < 1.7, 



10.9) 



with, as expected, Scenario 4 (green points) able to achieve these bounding values more 
easily than Scenarios 3 (red points) and 6 (brown points). Such departures from unity 
could easily be tested in view of a theoretically clean character of (3.20). 

Furthermore, in Fig. 19 we show the ratio of sin 2/3f over sin(2/3 + 2^ B J function 
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of 8f 3 for the scenarios considered. Similarly to r, the departure of this ratio from 
unity measures the violation of a golden MFV relation (3.11), this time between the 
CP-violating phases in the K — > nuu system and in the B® — B® mixing. As ipB d is 
constrained by the measured S^k s asymmetry to be at most a few degrees [40,68], large 
violations of the relation in question can only follow from the K — > -kvv decays. As seen 
in Fig. 19, they can be spectacular. 

10.12 The dependence on Sf 3 

As seen in Figs. 18 and 19 there are two oases in the values of 5f 3 

-10° < 5f 3 < 50°, 170° < 5f 3 < 250°, (10.10) 

with the desert between the two oases visibly but not densely populated. The origin of 
the oases is the constraint from S^k s - It is interesting to observe again a clear separation 
between Scenarios 4 and 6. While in Scenario 4 (green points) the first oasis is dominantly 
populated, in Scenario 6 (brown points) the second oasis is occupied. Moreover in the 
latter case there is an interesting excursion of points into the desert up to Sf 3 = 290°. 
The desert, instead, in which large effects simultaneously in K and B systems can be 
found, is dominated by red and blue points corresponding to Scenario 3 and the general 
scan, respectively. 

10.13 Comparison of Various Scenarios 

The plots in Figs. 10-19 are self-explanatory. Yet, we would like to make a few general 
observations: 

• There is a very clear distinction between scenarios 4 and 6. Scenario 4 can be 
considered as "£> s -scenario" as it gives interesting effects in the B s system. 

• On the other hand Scenario 6 can be considered as "if -scenario" as it admits 
spectacular effects in K decays. 

• Moreover, for certain sets of the LHT parameters obtained in the general scan, large 
departures from the SM predictions in K and B d s decays can be simultaneously 
found. 

10.14 B ttK Decays in the LHT Model 

We have finally investigated the impact of new physics contributions on B — > ixK decays 
that, for some time, signaled the presence of enhanced electroweak penguin contributions 
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with new large CP-violating phases. In a simple new physics scenario in which the 
universality between &s-penguins relevant for B — y ttK and srf-penguins relevant for 
K — y -nvv has been assumed [37], the enhanced electroweak penguins required to fit 
the B — y nK data implied spectacular effects in the K — » nuu system, similar to those 
shown in Figs. 10-12. 

Meanwhile, the data on B — > nK decays significantly changed [76-78] so that the 
SM predictions for the so-called R n and R c ratios [37] are nowadays within one standard 
deviations from the data. Consequently, within the new physics scenario of [37] no large 
effects in K — y nuu are expected. 

The picture is quite different in the LHT model, where the universality between bs 
and sd systems can be strongly broken. Following the approach of [37] to determine the 
hadronic parameters of the B — y nK system from the B — y 7T7T data and calculating the 
electroweak penguin contributions to B — y nK decays in the LHT model, we find that 
the new physics effects in these rare decays are smaller than the theoretical uncertainties 
present in non-leptonic decays. These small effects follow primarily from the smallness 
of complex phases in the bs penguins as given in (10.7) to be compared with ±90° taken 
in past B — y nK analyses. 

In summary, in the LHT model the smallness of new physics in B — >■ 7r.fr decays does 
not imply its smallness in K — y nuu decays as seen in Figs 10-14. 

10.15 What if the sin 2(3 problem disappears? 

Until now, our analysis was based on the tree level determination of the angle (3 that, due 
to the high value of Rb, is larger than the one measured through the S^k s asymmetry. 
It is conceivable that future tree level determinations of Rb will yield lower values for 
Rb and consequently for (5, so that the sin 2(5 problem will disappear. We have repeated 
the whole analysis for such a scenario, finding that the sin 2(3 problem has a very minor 
impact on our analysis, in particular in B s and K decays. For instance, large effects in 
K —y nuu decays and simultaneously in can still be found. 

10.16 What if the B d and B s Decays are SM-like? 

Finally, we can consider a very pessimistic scenario where BaBar, Belle and LHC will 
confirm all the SM expectations in B^ and B s decays, finding in particular both the S^k s 
and asymmetries very close to the SM values. This could already happen at the end 
of this decade. The question then arises whether in such a situation we could still expect 
large departures from the SM values in rare K decays, whose precise measurements will 
be available only in the next decade. It is evident from Figs. 10-14 that even in this, 
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pessimistic for .B-physics, scenario, spectacular effects in K — > txvv and also large effects 
in K L — > 7i°£ + £~ will be possible in particular in Scenario 6. 

11 Summary and Outlook 

In this paper we have analyzed for the first time rare K and B decays in the Littlest Higgs 
model with T-parity [I4-f6,29]. Together with our previous work [33] on observables 
related to particle-antiparticle mixing and B — > X SiC fy, the results of the present paper 
allow us to obtain a general description of FCNC processes in this model. 

On the technical side, we have presented a complete set of Feynman rules for the 
LHT model including also vertices with Goldstone bosons. The inclusion of 0(v 2 / f 2 ) 
corrections to some of the vertices has been performed here for the first time. These 
Feynman rules allow the calculations of 0{v 2 / f 2 ) contributions in arbitrary gauge and 
should turn out to be useful also for other observables. 

Using these rules we have calculated in the LHT model the short distance functions 
Xi, Yi and {i = K,d,s). In the LHT model these functions are complex quantities 
and carry the index i to signal the breakdown of the universality of FCNC processes, 
valid in the SM. The new weak phases in X h Y { and Z i: which are absent in the SM 
and models with MFV, imply potential new CP-violating effects beyond the SM ones. 
We would like to emphasize that the new parameterization of rare decays in terms of 
non-universal functions Xi, Yi and Zi can be applied to any model with new flavour and 
CP-violating interactions but the same low energy operators. 

With the functions X i: Y,- L and Z,- L at hand, one can straightforwardly calculate the 
branching ratios for a number of interesting rare decays. In particular, we analyzed 
K+ ->. vr+z/z/, K L -> 7T°z/z/, B s4 ->■ /2+ fi~ , B -> X s4 vv, B ->■ X s e + e~, K L ->■ n°£ + £- 
and B — > ttK. At all stages of our numerical analysis we took into account the existing 
constraints from electroweak precision studies [30] and from particle-antiparticle mixing 
and B — > X s ^ studied by us in [33] and from B — > X s i + i~ calculated here. 

The main messages of our paper are as follows: 

• The most evident departures from the SM predictions are found for CP-violating 
observables that are strongly suppressed within this model. These are the branch- 
ing ratio for K L — > n°vu and the CP-asymmetry S^. 

• Large departures from SM expectations are also possible for Br[K L — > tt°£ + £~) 
and Br(K + — >■ it + vu). 

• The branching ratios for B s4 — > and B — > X s ^vv, instead, are modified by 
at most 50% and 35%, respectively, and the effects of new electroweak penguins in 
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B — y irK are small, in agreement with the recent data. 

• Sizable departures from MFV relations between AM S:d and Br(B S:d — y fi + 'fi~) and 
between S^k s an d the K — y ixvv decay rates are possible. 

• The universality of new physics effects, characteristic for MFV models, can be 
largely broken, in particular between K and B s ^ systems. 

• The new physics effects in B — y X Sj( f) and B —y X s ^£ + £~ turn out to be below 5% 
and 15%, respectively so that agreement with the data can easily be obtained. 

One of the important findings of our paper is the presence of left-over singularities in 
the mirror sector that signals some sensitivity of the final results to the UV completion 
of the theory. This issue has been discussed in detail in the context of the LH model 
without T-parity in [28], is known from the study of electroweak constraints [29] and 
has been considered in Section 6 of the present paper. In estimating the contribution of 
these logarithmic singularities, we have assumed the UV completion of the theory not 
to have a complicate flavour pattern or at least that it has no impact below the cut-off. 
Clearly, this additional assumption lowers the predictive power of the theory. In spite of 
that, we believe that the general picture of FCNC processes presented here and in our 
previous paper is only insignificantly shadowed by this general property of non-linear 
sigma models. 

We conclude with probably one of the most important messages of this paper and of 
our previous analysis: 

• In spite of an impressive agreement of the SM with the available data, large de- 
partures from the SM expectations in B s decays are still possible. However, even 
if future Tevatron and LHC data would not see any significant new physics effect 
in these decays, this will not imply necessarily that new physics is not visible in 
Kl — y 7r°z/z/, K + — y -n^vv and — y n £ + £~. On the contrary, as seen in the 
case of Scenario 6, large departures in these three decays will still be possible. It 
may then be that in the end, it will be K physics and not B physics that will offer 
the best information about the new phenomena at very short distance scales, in 
accordance with the arguments in [79,80]. 

Acknowledgements 

We would like to thank Bjorn Duling for checking the Feynman rules given in Appendix 
B and William A. Bardeen, Ulrich Haisch and Kazuhiro Tobe for useful discussions. This 
research was partially supported by the German 'Bundesministerium fur Bildung und 
Forschung' under contracts 05HT4WOA/3, 05HT6WOA. 



48 



A Non-leading contributions of <J> 



In this Appendix we show that the scalar triplet $ does not contribute to the decays 
considered in the present paper at 0(v 2 / f 2 ). 

In principle, also the scalar triplet $ could contribute to the decays considered in the 
present analysis. The corresponding diagrams can be obtained from the ones shown in 
Figs. 1 and 2 by simply replacing W% by ± and Z H , A H by 0°, P . Therefore we also 
derived the Feynman rules for the vertices containing $. We found them to have the 
following generic structure: 

where c, d are 0(1) coefficients depending on which qjj, q and component of $ are 
involved, and Pr,l — (1 ± 7s)/2 are the usual projectors, qu and q denote mirror and 
SM fermions, respectively, and Gr, Gl are heavy and light gauge bosons. 

As we have set the masses of the external quarks and leptons to zero throughout our 
analysis, we obtain 

q H $q ~ O (^p) and $G H G L ~ O (^jj . (A.2) 

It can now be easily seen that each diagram contains at least two such vertices, so 
that they are suppressed by 0(v 4 / f A ). 

Furthermore, in contrast to the LH model without T-parity, all particles running in 
the loops have O(f) masses, so that no cancellation of v/f factors due to large mass 
ratios, as encountered in [26,28] for diagrams with T + exchanges, can appear. 

Thus we find that $ contributes to the decays in question - and more in general to 
all decays with external SM fermions - only at 0(v 4 /f A ). For this reason we do not give 
any Feynman rule for the interactions of the scalar triplet $. 



B Relevant Feynman Rules 

In this Appendix we list all Feynman rules relevant for the analysis performed in the 
present paper, and describe briefly how they have been derived. We note that given 
the Feynman rule for a vertex, the rule for the conjugate one can be obtained through 
the replacement (vertex)— > —(vertex)'''. There follow, in particular, the prescriptions 
Pl — > —Pr, Pr — > —Pl and 7 m Pl,/? — > 1^Pl,r- A similar, but more detailed description 
for the LH model without T-parity can be found in [17] and in [24,28], where some of 
the Feynman rules given in [17] have been corrected. 
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B.l Fermion— Gauge Boson Couplings 

The fermion-gauge boson couplings can easily be obtained from the fermion kinetic 
term [41] 

fermion = ^l^^l + *2^^*2 + ^ D% + ^ D% , (B.l) 

where 

D\ = d, - VligQWl ~ sftig'Y^B^ - V2ig'Y^B 2 , , (B.2) 
Dl = d, + V2igQfW^ - V^ig'Y^B^ - yfagfY™ , (B.3) 
D'; = d,- V2ig'Y^B^ - V2ig'Y^B 2 , , (B.4) 

by inserting the mass eigenstates of all particles involved. The U(l) charges can be 
obtained from gauge invariance of the Yukawa couplings and T-parity [41] and are given 
in Table 2. 
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Table 2: U{\)\ x U{1) 2 quantum numbers of the fermion fields. 

Similar kinetic terms can be written down for the right-handed SM fields and for 
t'iRi i = 1? 2. 

The kinetic term for the right-handed mirror fermions can be written down by using 
the CCWZ approach [31,81] and is given in explicit terms in [41]. As the mirror fermions 
are purely vector-like under SU(2)l x U(1)y, their couplings to SM gauge bosons are 
proportional to 7 M . Thus there is no need to consider the CCWZ kinetic term any 
further, as the couplings of the right-handed mirror fermions to SM gauge bosons have 
to be equal to the ones of the left-handed mirror fermions. 



Fermion couplings to SM gauge bosons 
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The couplings to the heavy gauge bosons have already been given in [32] . We confirm 
the findings of these authors, and include also v 2 /f 2 corrections to the couplings of 
the neutral gauge bosons. These corrections turn out to be irrelevant in the decays 
considered in the present paper, but could be of relevance for other processes. While 
our paper was being completed, the rules for the quark sector appeared also in [42], but 
without performing the v/f— expansion and without taking into account the flavour 
mixing matrices Vh u and Vnd- 

Note that the couplings to the heavy gauge bosons are purely left-handed as there 
is no kinetic term including both heavy and light right-handed fermions. The only 
exception are the couplings involving T_ together with T + or t. 
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gauge bosons 


u l H A H u j 


(_W _ i_l r „^ 
\ 10 2 XH p , 


) (v H uhrp L 


u % H Z H u 3 


(ig _ igL TlJ ^l\ 


(VHuhrPL 


d L H A u d? 


(—iaL + i9 r ^ 


) (Vh^Pl 




( ig ig' ry-r T V2 ' 

\ 2 10 XH p~ y 


) (VndhrPL 



51 







u % H Z H t 


m _ ( ij_ r 2 1 i£ T[J 

2 1 4 J/ i t 10 H 


(v Hu ) i3 rp L 


u l H A H l + 




u l H Z H T + 


i {x L j(v Hu ) i3 rp L 


f_A H t 


-\ig'Y (x L jP L + JxZ (l - 


-${l-x L ){\-x L ))P R ) 


f-Z H t 




T_A H T + 




xl{1 + $x l (1-x l ))Pr 


f_Z H T + 


iig'xH jii^PL + Vi - x l p r ) 




js(v Hd ) 


■ -V* Pt 


ti Jri 






11 11 








^(V Hu hx L f^P L 


v\jA H vi 


( iff' iff v 2 \ 
I 10 ~~ ~2 X HJ2 j 




uirZuU'j 


( id _i_ iff' ™ T v 2 \ 
\ 2 ^ 10 XH J^ J 




t H A H V 


( ig' J_ ig r,, v 2 } 






\~2 + W XH J 2 


) (V^-t^l 




^(Vh^Pl 









B.2 Fermion— Goldstone Boson Couplings 

The fermion couplings to T-even and T-odd Goldstone bosons can be derived from the 
Yukawa interactions (2.23), (2.38) and (2.40) and from the mass term (2.20). Note that 
(2.40) with X = (S 33 ) _1 / 4 corresponds to case A in [43]. Using instead X = (S33) 1 / 4 
(case B in [43]) does not modify the fermion couplings to Goldstone bosons as required 
by gauge invariance. 

As these terms have again to be expressed in mass eigenstates of the particles involved, 
we have to take into account the 0(v 2 /f 2 ) mixing of the Goldstone bosons and scalars 
[30]. Following the steps of [30], we find the mass eigenstates, denoted here with a 
subscript "to", to be 

< = *» (1 - 4p) • < R5 > 
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h m = h, (B.7) 

<t>° m = 0° (l - ^) , (B.8) 

C = <P P ~ (VWr, ~ V2co° + P ) ^ , (B.9) 

^ = ^{}-^p) tiu± W 2 ' (R10) 

0+ + = ++ , (B.ll) 

^ = r? 24? (R12) 

o o %(^° + 4y^0 p ) - 4VE V (5g + 6<fo g ) t; 2 

<"m = <" ^ ? , (B.13) 

^ = w± ( 1 -2&)^H^f (R14) 



In what follows, we will drop the subscript "to" , as it is clear that all Feynman rules are 
given in terms of mass eigenstates. 

Note that the mixing of Goldstone bosons and scalars in question affect the corre- 
sponding Feynman rules at 0(v 2 /f 2 ), thus without considering it, one would not be 
able to get the right 0(v 2 /f 2 ) corrections. Also the 0(v 2 /f 2 ) corrections to the particle 
masses have to be taken into account. 

We would like to caution the reader that this mixing has not been taken into account 
in the derivation of the Feynman rules for $ given in [41]. 



Fermion couplings to SM Goldstone bosons 




V2M Wl (< Pl m d P R) (Vckm)^ 


in + d j 




f + 7T + d j 


Jm Wi , 7 i 1 + P d2 ) PL < P *_ ^KM) t , 




9 m u „ X 
2Mz L cose w '5"ij 


tn°t 


9™i ( 1 r 2 f 2 \ 
2M Zl cos6» w ^ X L p ) 15 


T + TT°T + 


9' m T + J2 v 2 „, 
2M Zl cos 9 W X L pl5 


f + n°t 


iMz^J m T+ (l + f(d 2 - x i))p L -m t P R 


d l ^d? 


9m%d ^S- 

2M Zt cos6 w > 5Ul 3 
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ui I 7T + d 3 H 


9 v2 m d fr 
8V2M Wr 1 J 








9ml 

2M Zr cos 9 W ' 5 "*J 


£W 


2M Zr cos9 w l 5U V 




9 v 2 e X 
8V2M Wr P m Hi°i] 



The leading order contributions to the couplings of fermions to the heavy Goldstone 
bosons have already been given in [32]. We included also 0{y 2 / f 2 ) corrections, necessary 
for the calculation of rare decays, and the contributions proportional to the SM fermion 
masses. We note that our Feynman rules for these vertices differ by sign from the ones 
given in [32], as our sign conventions for the mass and Yukawa terms (2.20), (2.23), 
(2.38) and (2.40) are chosen such that we obtain positive and real- valued masses, as also 
done in [43]. 



Fermion couplings to heavy Goldstone bosons 


u l H u + d j 


V2M Wff ( m Hi P L <P 


r) (VHd)ij 


u l H u°u j 


9 

2M Z[I 


m%{l + $(\ ta X)) 


Pl ~ m{P R 


( V Hu)ij 




9 

2M Zh 


u / 1 , v 2 ( 1 x„ x 2 L 
m Hi [ L + p [ 8 t<m0 w 2 


))P L -m t P R 


(VHu)i3 




iff^f {m u m P L ~ m T+ P R ) (V Hu ) i3 




gm t v p 
2M Zu r R 




gm T+ v 2 / ixH 
2Mz H f 2 \ l tandw 






9' 

10M Aff 


m m f 1 + 75 (| + Zfz-tanflw 


) P L - m{P B 


(V Hu )ij 




9' 

iom Ah 


m Hi U + jz (| + x H taja.9 w - 


x -L^P L -m t P R (V Hu ) i3 






Pr) (Vnu)i3 


f_T]t 


2g'm t f ( v 2 ( x'l , 1 \ \ p 
5M A „ v\ L p \ 2 + 6j) r R 




2g'x L m T+ 1 v 2 , 3 , 
5M Aff \ L p \2 X L 


2x l + 1))Pr 


d l H oj~ui 


V2M W „ m »A l W)^ 


m{P R {V Hu )ij 


d l H uj~t 


V2M Wh m Hi ( X p ( 8 + 2 )) " 


Pl ~ m t P R 




d i H oj-T + 


V2M Wff XL J «* Pl ^ 


^Pr) (Vnu)i3 


d^M 


9 

2M Zn 


m m (l + p ( 4 + tan0 w ) 


) P L - m? d P R 


{V H d)ij 
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9' 

iom Ah 


m d Hi (l-$(l + x H tfm8 w )] 








V2M W „ ™M PHI)* 




9 

2M Zr 


m Hi( 1 + j^( 8 tan^)) 


P L -miP R {V Hv )ij 




9' 

10M Aff 


m m + (l + x H tfw6w) 


) 


P L - miP R 






V2M W „ m Hi i 1 8p) P L 


- m 




t H u°e j 


9 

2M Zr 


m Hi (i + p ( 4 + tan^, 





P L -m{P R {V H i)n 


t H T]l j 


9' 

iom Ar 




: 







B.3 Triple Gauge Boson Couplings 



The self-interactions of the gauge bosons can be obtained, as usual, from the gauge boson 
kinetic term 



-gauge 



-^Tr(F^P^) - iTr(P^P^) - \b\ v B^ - \b%B^ (B.15) 



by inserting the mass eigenstates of the gauge boson fields. Here, F l and B l denote 
the field strength tensors of the gauge groups SU(2)i and U(l)i (i = 1, 2), respectively. 
The momenta k, p and q are defined to be incoming. 



Gauge boson self-interactions 


W^{k)W^{p)Z p L {q) 


ig cos 9 W [g pu (k - p) p + g up (p - q) p + g p »(q - k) u ] 


W^{k)W^{p)A p L {q) 


ie [g pu {k - p) p + g up {p - q) p + g p ^(q - k) u ] 


W^(k)Wr( P )Z p H (q) 


ig [g^(k - p) p + g up (p - q) p + g p ^{q - k) u ] 


W^(k)W^(p)K(q) 


igxujz [g^ u (k - p) p + g up (p - qY + g pp (q - k) u ] 


W^(k)W-(p)Zi(q) 


ig cos 9 W [g pu {k - p) p + g up {p - q) p + g pp {q - k) v ] 


W^{k)W H ^p)A p L {q) 


ie \g pv {k - p) p + g vp {p - q) p + g m {q - k) u ] 



B.4 Triple Gauge Boson— Goldstone Boson Couplings 

The kinetic term for the non-linear sigma model field S is given by 

C = ^Tr [(D^ (D"E)] , 
where the covariant derivative is defined through 



PVE = d/jT, -V2iJ2 Wl (QJE + SQf ) + g'B m ft-E + Elf)] . 

3=1 



(B.16) 



(B.17) 
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From this term the interactions of the Goldstone boson fields with the SM and heavy 
gauge bosons can be obtained. The mixing of the Goldstone boson and scalar fields, as 
described in Appendix B.2, has also to be taken into account. 
All momenta are again defined to be incoming. 



Gauge boson-Goldstone boson interactions 








/ 9 \ 

-gM ZH (l-^)g^ 


W^Z v l uj- 


9 m Wh cos e w [i 4/2 ; os2 ,J<r 
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-eM WH g^ 
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-eM WL g^ 


W^A h tt- 
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ie (p - qY 
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~i 9 i (P - if 


uj + (p)io°(q)W- 11 


- ig (i-^yp- q y 


u + (p)r](q)W^ 


1 2iVl f 2[ P q > 


^(p)co%q)W^ 




n+(p) V (q)W H » 


-i^ f ip-qY 


n+(p)n°(q)W^ 


-il{p-qY 



C Relevant Functions 



In this Appendix we list the functions that entered the present study of rare and CP- 
violating K and B decays. Both the SM contributions and the new physics contributions 
coming from the T-even and T-odd sectors are collected. The variables are defined as 
follows: 



m 2 Wl 



Tfln 



X T 



(q = c,t), 



m 



Hi 



m 



Hi 



Ah 



Zi a with a 



tan 2 9 W ' 



(2 = 1,2,3), 



(C.l) 
(C.2) 
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m 



y = 

-^SM(^t) 

Y SM (x t ) 

ZsM(x t ) 



Hi 



rn 



HI 



M 2 Wh MI h ' 



y =ya, i] = 



8 
8 



x t + 2 3x t 
x t 

x t 



6 



1 ' (a^-1) 2 
4 3x t 



log x t 
log a; t 



x t - 1 (x t - l) 2 

18x4 - 163x3 + 259x 2 - 108x t 



--logXi + 



32x 4 - 38x 3 - 15x 2 + 18x t 



CM 
D (y) 
E (y) 
D' (y) 
K(y) 



72(x t - l) 4 



y-6 3y + 2 
y-l (y-1) 2 



lU(x t - If 
logxt- 



4 -19y 3 + 25y 2 y 2 ^/ 2 - 2y-6) 

-9 hgy+ 36(y-l)3 + 18(y-l) 4 l0g2/ ' 

2 j/ 2 (15-16g/ + 4y 2 ) l y(18-lly-y 2 ) 
"3 l0g2/ + %^l7 ^ + 12(1 -y) 3 

(3y 3 - 2y 2 ) (8y 3 + 5y 2 - 7y) 

logy + 



2(y-l) 4 
3y 2 



2(2/- 1)^ 



logy + 



12(y - 1)3 

(y 3 - 5y 2 - 2y) 
4(y - l) 3 ' 



(C.3) 

(C.4) 
(C.5) 



-3 + 2x t -2x 2 x t (-A-x t + 2xf)\ogx t , (3 + 2x t )logx T 

XT) = g( -i + x,) sFTT^F + 8 

(3 + 2x t - 2x 2 ) x t (2-x t + 2x 2 )logx t (3 + 2x t ) logx T 



Va(a; t ,xr) 



8(-l + x t ) 



8(-l + x t ) 5 



(C.6) 

(C.7) 
(C.8) 
(C.9) 
(CIO) 
(C.ll) 

(C.12) 
(C.13) 



(D')sm = D' (x t ) , 
(£>')lht 



(D )even + t s ^T D i , 



It) 2 



T-iodd 



#odd(*) 



4/ 2 



&\D' odd (z2) -D' odd (zi)) +es , (D' odd (z 3 )-D' oM (z 1 )) 



( s ) f n' 



(C.14) 
(C.15) 

(C.16) 
(C.17) 
(C.18) 
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R2(Zi) = 

F 2 (z i ) = 



Zj log Zj _L 

(1 - z,f + 1- Zl 

zf log Zi 1 



(i - Zl y i - zi 



(C.19) 
(C.20) 



F vV { Zi ,y-W H ) = -z i -F 5 (z i ,y)-7F 6 (z i ,y)-9U(z i ,y) , 
F^( Zl ,y;W H ) = ^z i -F 5 (z i ,y)-7F 6 (z i ,y) + 3U(z i ,y) . 



(C.21) 
(C.22) 



F 5 (zi,y) 
F e (zi,y) 
U (zi,y) 



zf log Zi 
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y 3 logy 



- Zi) (y - z^ {l-y){zi-y)' 
zf log Zi y 2 log y 



+ 



(1 - z^ (y - z^ (1 - y) (zi - y) 
zf log Zi y 2 log |/ 



+ 



(zi - y) (1 - *) 2 (v - *) (1 - l/) 2 (! - *) (1 - V) 



(C.23) 
(C.24) 
(C.25) 



G(zi,y;Z H ) = --U (z,y) , 
G l (zi,y>;A H ) = ^G(z' l ,y' l ;Z H ) , 
G 2 (zi,y;r)) 



zf log Zi 



+ 



10a (1 - z^ (rj - z^ (zi - y) 
y 2 log y 



+ 



?/ 2 log f] 



(1 -y)(rj- y) (y - zt) (1 - 77) (z { - rj) (77 - y) 
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